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We present the first in a series of papers studying with Chandra the X-ray properties of a sam- 
ple of 28 early-type galaxies which span ~3 orders of magnitude in X-ray luminosity (L x ). We 
report emission-weighted Fe abundance (Zp e ) constraints and, for many of the galaxies, abundance 
constraints for key elements such as 0, Ne, Mg, Si, S and Ni. We find no evidence of the very sub- 
solar Zp c historically reported, confirming a trend in recent X-ray observations of bright galaxies and 
groups, nor do we find any correlation between Zp c and luminosity. Except in one case we do not find 
evidence for a multi-phase interstellar medium (ISM), indicating that multi-temperature fits required 
in pr evious AS CA analysis arose due to the strong temperature gradients which we are able to resolve 
with Chandra\ We compare the stellar Zp e , estimated from simple stellar population model fits, to 
that of the hot gas. Excepting one possible outlier we find no evidence that the gas is substantially 
more metal-poor than the stars and, in a few systems, Zp e is higher in the ISM. In general, however, 
the two components exhibit similar mctallicitics, which is inconsistent with both galactic wind models 
and recent hierarchical chemical enrichment simulations. Adopting standard SNIa and SNII metal 
yields our abundance ratio constraints imply 66±11% of the Fe within the ISM was produced in SNIa, 
which is remarkably similar to the Solar neighbourhood, and implies similar enrichment histories for 
the cold ISM in a spiral and the hot ISM in elliptical galaxies. Although these values are sensitive to 
the considerable systematic uncertainty in the supernova yields, they are also in very good agreement 
with observations of more massive systems. These results indicate a remarkable degree of homology 
in the enrichment process operating from cluster scales to low-to-intermediate Lx galaxies. In addi- 
tion the data uniformly exhibit the low Zo/ZMg abundance ratios which have been reported in the 
centres of clusters, groups and some galaxies. This is inconsistent with the standard calculations of 
metal production in SNII and may indicate an additional source of a-element enrichment, such as 
Population III hypernovae. 



Subject headings: Xrays: galaxies- 
galaxies: halos- 

1. INTRODUCTION 



- galaxies: elliptical and lenticular, cD- 
galaxies: ISM 



galaxies: abundances- 



The entire history of star-formation and evolution 
leaves its chemical signature in the hot gas of early- 
type galaxies. X-ray observations therefore provide a 
natural and po werful diagnostic tool to unlock this in- 



formation (e . g. Loewenstein fc Mathews) 



Brighcnt| |2003| 
surements of intersl 



However, historica 



1991 



Mathews 



ly X-ray mea- 



stellar medium (ISM) abundances have 
been, probl ematical, ao typified by the go called "Fe dio 
crepiacj 



ray observations typically find Zp c ~0.3-0.5. The metal 
content of the intra-cluster medium (ICM) is primarily 
attributed to the stellar ejecta from giant elliptical galax- 
ies, and so it is difficult to envisage there being lower 
metal abundan ces in individual galaxies than in the ICM 
flRenzini|[i99"7l) . 

Attempts to reproduce the low Zp e in early-type galax- 
ies from a modelling standpoint by, for example , incorpo- 
rating ongoing ac cretion of unenriched gas (e.g. Brighcnti 



observations of these galaxies tended to imply extremely 
sub-solar metal abundances (generally expressed as Zp c , 
the Fe abundance with respect to the adopted solar stan- 
dard, since Fe has th e strongest diagnostic lines in th e 
soft X-ray band) (e.g. Loewenstein & Mushotzky 1998), 
in stark contrast to ~solar abundances in the stellar pop- 
ulation. Since individual galaxies are not "closed boxes" , 
the ISM is believed to be built up primarily through 
stellar mass-loss and Type la supernovae (SNe) ejecta. 
These two crucial ingredients lead classical "wind mod- 
els" of gas enrichment to predict highly super-solar Zpa 



omo.tical. au tvpmoa by tnc oo called fc dip n , m l 71 ,. , — y 1 p T — 

QArimoto et al |h997p Early Rasat and ASCA fc Mathewsj 1999|) or allowing complex star-formation his 



tories (e.g. Kawata & Gibson 2003 ) have not been en- 
tirely successful, le ading some to point a finger at the 



spectral- modelling. Arimoto et al. ( 1997 ), for instance, 



called into question the plasma codes being fitted to the 
data, particularly in light of uncertainties associated with 
the Fe L-shcll transitions in X-ray emiss ion plasma codes 
This p oint was investigated further by Matsushita et al 
( 2000 [), who argued that, having degraded the data qual- 



in these galaxi es (e.g. Ciotti et al. 1991; Loewenstein & 
Mathews! |1991| ). The problems caused by this discrep- 



ity in the vicinity of the Fe L-shcll, Zp c > 0.5 can be found 
in the highest X-ray luminosity (Lx) galaxies. Consistent 
results from spectral-fitting with plasma codes that treat 
the Fe L-shell transitions differently would seem, how- 



ancy are exacerbated further when attempting to under- 
stand gas enrichment in clusters of galaxies, for which X- 

1 Department of Physics and Astronomy, University of California 
at Irvine, 4129 Frederick Reines Hall, Irvine, CA 92697-4575 



ever, t o conflict with this explanation (e.g. Buote et al 
2003b| ). Although remaining errors in the treatment of 



the Fe L-shell in the plasma code s can be important for 
high-resolution spectroscopy (e.g. Behar et al. 2001), at 
CCD resolution the effects are substantially washed out 
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Humphrey & Buote. 



so that they contribute only a ~10-20% sy stematic un 



certainty to the abundance measurement (Buote et al 

Eoo3tfL 



Pe 'haps a more natural alternative was suggested by 



Buote & Fabian (1998), who demonstrated that fitting 
a single temperature model to the emission spectrum of 
intrinsically non-isothermal hot gas gives rise to a sub- 
stantially under-estimated abundance (not to mention a 
significa ntly poorer fit ), an effect dubbed the "Fe bias" 
(see also Buote _ 2000b|). This effect had previ ously been 



recogn ized by Buote fe Canizarcs ( f 994 ), and Trinchicri 
et al] ( 1994 ) found that fitting a two-temperature model 
to the Rosat spectrum of the bright elliptical NGC 4636 
resulted in poorly-constrained abundances which could 
be consistent with solar values, in contrast to low Zp 
required for single-temperature fits. A spectrally hard 
(kT> 5 keV) component had long been recognized in the 
ASCA spectra of the many early-type galaxies, which 
was attributed to emission f rom undetected X-ray bina- 
ries ( Matsushita et al. 1994 ). The composite spectrum 
from these sources, which dominate the emission of low- 
Lx galaxies, can typically be approximated as a pure 
bremsstrahlung model. Therefore, fitting only a single- 
temperature hot gas model to the spectra of such galax- 
ies a lso tends to give rise to unphvsically l ow abundances 



(e.g. Fabbiano et al. 1994 ; Kim et al. 1996 , who found es- 



sentially unconstrained abundances when incorporating 
a term to account for this effect in their Rosat analy- 
sis). This effect is clearly disti nct from, albei t related to, 
the Fe bias (see discussion in Buote 2000b ). Although 
the spectral-shape of the unresolved source component 
is sufficiently hard that it cannot move to mitigate, even 
in part, the Fe b ias arising from multi-t emperature hot 
gas components, Buote fe Fabian ( 1998 ) still found Zp c 
consistent with solar in the lower-Lx systems, in which 
the fit implied only a single hot gas component, plus un- 
detected sources. However, the poor signal-to-noise ra- 
tio (S/N) characteristic of the lower-Lx galaxies tended 
to produce poorly-constrained abundances, so that very 
sub-solar abundanc es could not be ruled out. 

Recent Chandra and XMM measurements of abun- 
dances in X-ray bright galaxies and the centres of groups 
have tended to support ^solar or even slightly super - 



[Tamura et al. 


2002; 


Gastaldello & Molcnd: 


^002 


; Buote 


P002L Ku et al 


2002 


'Sullivan et al. 2003]; |Kim & Fab- 



bian| 2004b| ). In contrast, very sub-solar Zp ar e still 



being repo r ted in the lowest Lx /Lr systems (e.g. Irwin 
et al] 2002 ; Sarazin et al. [2001 ). Perhaps the most dra- 
matic examples of this la tter effect are the three very 
low-Lx galaxies for which O'Sullivan & Ponman (2004) 
reported Zp c <0.1. This apparent lack of consistency m 
the enrichment processes operating on different scales is 
intriguing, and it remains to be assessed whether it is an 
artefact of the poorer S/N or a real effect in low-Lx sys- 
tems. Although a major problem for our understanding 
of galaxy evolution, a full understanding of this effect 
is inhibited by the lack of interesting abundance con- 
straints in galaxies with Lx intermediate between these 



two extremes. In a recent paper, Humphrey et al. (2004), 
we made some initial progress in this area by report- 
ing constraints for the normal, moderate-Lx SO galaxy 
NGC 1332 and the elliptical NGC 720, in both cases 
strongly excluding the extremely sub-solar (Zp < 0.4) 



abundances historically reported for these systems. Cou- 
pled with similar constraints on the abun dances in the 



moderate L x/Lr radio galaxy NGC 1316 ( |Kim fe Fab 



biano 2003), this would hint at a consistent picture of 
enrichment from cluster to moderate-Lx galaxy scales. 

In addition to the insight into enrichment afforded from 
global Zp e measurements, the a-element abundance ra- 
tios, with respect to Fe, provide an additional power- 
ful diagnostic. Different types of SNe inject material 
imprinted with characteristic chemical "fingerprints" , so 
that the measured abundance pattern can be used to as- 
sess the relative contribution of SNIa and SNII to ISM 
enrichment. Early attempts to extract this information 



(e.g. Matsushita et al 



were largely affected by a 



2000 



Finogucnov & Jones 200C) 
ailure to treat the Fe bias, 



although a significant contribution of SNIa to the en- 
richment process was indica ted. Mor e recently, studies 
with high-quality XMM and Chandra data have revealed 
typically ~70-90% of the ISM enrichment in the cen- 
tres of groups and clusters seems to have its origin in 
SNIa, which is close to the ^75% value in the Solar neigh 



bourh ood (e.g. |Gastaldcllo fe Molendi| |2002| ; [Buote et al 
2003b| ). For NGC 1332 and NGC 720, we were able to 
obtain interesting constraints on the a-to-Fe ratios in 
some of the lowest-Lx systems to date, again inferring 
~ 70-80% of the Fe to have been produced in SNIa. This 
also supports the suggestion of homology in the enrich- 
ment process over different mass-scales, at least down to 
moderate-Lx galaxies. 

Although these results are intriguing, it is by no means 
clear that the two galaxies we considered are representa- 
tive, nor is it clear that all X-ray brig ht galaxies follow 
the t rend of ^solar abundances (e.g. Sambruna et al. 
2004 ). Furthermore, these studies have not elucidated 
the processes which may be giving rise to very low-Zp c 
measurements in the faintest systems. In light of these 
results, therefore, the next logical step is to consider a 
uniformly- analyzed sample of galaxies spanning a large 
range of Lx, and particularly expanding the number of 
moderate-Lx galaxies studied. In this paper, we present 
the results of a study of metal enrichm ent in a sa mple of 
28 early-type galaxies drawn from the Chandra archive. 
In subsequent papers we will discuss the gravitating mass 
and point-source populations of these objects as a whole. 
The galaxies have been carefully selected to span the 
available X-ray luminosity range, from group-dominant 
to low-Lx galaxies. 



In many respects, Chandra ACIS is the natural instru- 
ment with which to undertake such a study. Although 
XMM has a significantly higher collecting area, the anal- 
ysis of faint, diffuse sources i s complic ated by difficulties 



in treating the background. Chandra also has an intrin- 



sic advantage in being able to resolve out a substantial 
fraction of the X-ray binary contribution into individual 
sources, which must otherwise be d isentangl ed spectrally. 
The excellent spatial resolution of Chandra also provides 
an unprecedented opportunity to investigate any spatial 
temperature variation, which would provide the natural 
source of the Fe bias. Imaging spectroscopy at CCD reso- 
lution is well-suited to determine reliable abundances in a 
gal axy. In fa ct there arc a number of drawbacks to using 
the Chandra and XMM gratings instead for such a study. 
The extended nature of the sources makes grating spec- 
troscopy extremely challenging. The significantly lower 



ISM abundances in early-type galaxies. 
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effective area of the grating spectrographs in comparison 
to the ACIS CCDs and, especially for the XMM RGS in- 
strument, the more limited bandwidth both exacerbate 
this problem. Although one of the key issues of inter- 
est is the determination of multiple temperature compo- 
nents in the extraction aperture, this tends to produce 
features broad enough to be evident in CCD spectra. In 
fact, based on high S/N data of the group NGC5044 
there was excellent agreement in the abundances deter 
mined with the XMM gratings and the XMM and \Chan 



dra@CDS ( |Buote ct al.||2003b| ; fTarnura et alJ[200^ )7coE 



firming that CCD resolution is sufficient to obtain reli- 
able abundances. 

Throughout this paper we adopt the l atest sola r abun- 
dances standard of Asplund et al. ( 2004 , see § |4.l| ), which 
resolve several previously-noted discrepancies between 
"photospheric" and "meteoritic" values. All error-bars 
quoted refer to the 90% confidence region, unless other- 
wise stated. 

2. TARGET SELECTION 

By its nature, the Chandra\ archive contains a nonuni- 
form sample of galaxies. Although this prevents our 
choosing a statistically complete sample, we selected 28 
galaxies approximately spanning the range of measured 
Lx, from ~ 8 x 10 39 -10 43 erg s -1 . We only considered 
relatively nearby galaxies which had been observed for 
a total of at least 10 ks with the ACIS-S or ACIS-I, 
and without any grating employed. For selection pur- 
po ses X-ray luminosi ti es we re taken from the catalogue 



of p 'Sullivan et al 



(2001), although all luminosities 
ha ye s ubsequently been recomputed in the present work 
(§3^2). To an initial sample of 26 galaxies chosen in this 
way, we also added two further interesting objects not 
in the O 'Sullivan catalogue — NG C 1132. one of the clos 



est ex amples of a "fossil group" ( Mulchaey fc Zabludofi 



1999| , F. Gastaldello et al, 2005, in prep.) and the rela- 



tively isolated, moderate-Lx galaxy NGC 1700 which has 



an unusually h igh X-ray ellipticity, which Statler & Mc- 
Namara fl2002|) argued may indicate rotational support 



Two of our sample, N GC 1332 and NGC 720 have al 
ready been discussed in Humphrey et al. (2004), and we 
adopt the results from that work here. The sample con- 
tains 10 high-Lx (logioLx^41.5) galaxies, 12 moderate- 
Lx(logioLx— 40.5-41.5) and 6 low-L x (logi L x ^ 40.5 ) 
galaxies. Most of the sample was observed with ACIS-S, 
although in a few cases ACIS-I was employed. To give ex- 
tra coverage at large radii both ACIS-S and ACIS-I data 
for the two bright galaxies NGC 1399 and NGC 4472, 
were analysed together. 

A s ummary o f the properties of the galaxies and details 
of the Chandra exposures are given in Table [l]. In order 
to provide accurate luminosity estimates, we searched 
the literature for reliable distance estimates. Where pos- 
sible, we adopted those de termined from surf ace bright- 
ness fluctuations (SBF) by Fonry et alj (2001, corr ectin, 
for an improved Cepheid zero-point: jJcnscn et al 



2003") 



or Jensen et al. ( 2003 ). Alternatively, we used distance s 



determined from the D n — a relation ( Faber et al. 1989 ) 
or the redshift, corrected for Virgo-centric flow, as given 
in LED A. We assumed H =70 km s -1 Mpc" 1 . 

Of this sample, 6 systems (IC4296, NGC 507, 
NGC 741, NGC 1399, NGC 4472, NGC 7619) appear to 
be the central galaxies in substantial, optically-identified 



groups so their X-ray emission may be to some extent in- 
tertwined with that of a surrounding intra-group medium 
(IGM). In the present context, it suffices to consider that 
all the systems comprise a continuum over a range of 
mass-scales. In a subsequent paper, we will discuss the is- 
sue of group membership, and the total gravitating mass, 
in detail. 

3. DATA REDUCTION 



CI AC 3.1 and Hea- 



Chandra 



Fo r data reduction we used the 
sofjd. 'S software suites, in conjunction with 
CaUU\ calib ration database 2.28. For spectral-fitting we 
used |Xspec 11.3.1. In order to ensure the most up-to- 
date calibration, all data were reprocessed from the "level 
1" events files, following the standard Chandra data- 
reduction threads 2 . We applied corrections to take ac- 
count of a time-dependent drift in the satellite gain and, 
for ACIS-I observations, the effects of "charg e trans fer in- 
efficiency" , as implemented in the standard CI AC tools. 

To identify periods of enhanced background ("flar- 
ing" ) , which seriously degrades the signal-to- nois e (S/N) 



and compl icates background subtraction (e.g. Marke- 
yitch||2002D , we accumulated background lightcurves for 



each exposure from low surface-brightness regions of the 
active chips. We excluded obvious diffuse emission and 
data in the vicinity of any detected point-sources (see 
below). Periods of flaring were identified by eye and ex- 
cised. Any residual flaring which is not removed by this 
procedure will be sufficiently mild to have negligible im- 
pact in the centres of bright galaxies. However, in the 
fainter systems even very mild background variation can 
have a significant impact on our results, and so we treat 
these systems with extra care (§ |5.2j ). The final exposure 
times are listed in Table |l[ 



Point source d etection was performe d using the CI AC 
tool wavdetect ( Freeman et al. 2002 ). In order to im- 
prove the likelihood of identifying sources with peculiarly 
hard or soft spectra, full-resolution images were created 
of the region of the ACIS focal-plane containing the S3 
chip in the energy-band 0.1-10.0 keV and, so as to iden- 
tify any unusually soft or hard sources, also in the bands 
0.1-3.0 keV and 3.0-10.0 keV. Sources were detected sep- 
arately in each image. In order to minimize spurious 
detections at node or chip boundaries we supplied the 
detection algorithm with exposure-maps generated at en- 
ergies 1.7 keV, 1.0 keV and 7 keV respectively (although 
the precise energies chosen made little difference to the 
results). The detection algorithm searched for structure 
over pixel-scales of 1, 2, 4, 8 and 16 pixels, and the de- 
tection threshold was set to ~ 10 -7 spurious sources per 
pixel (corresponding to ~0.1 spurious detections per im- 
age). The source- lists obtained within each energy-band 
were combined and duplicated sources removed, and the 
final list was checked by visual inspection of the images. 
A full discussion of the point source populations will be 
given in a subsequent paper. In the present work, the 
data in the vicinity of any detected point source were 
removed so as not to contaminate the diff use emis sion. 
As discussed in [Humphrey fc Buote ( 2004 , see also 



P 



Fabbianc 



2004. 



Kim 



al) a significant fraction of faint X-ray 
II not have been detected by this pro- 



mary sources wil. 
cedure, and so we include an additional component to 



http: / /cxc. harvard.edu/ciao/thrcads/indcx. html 
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TABLE 1 

The sample 



Galaxy 


Type 


Dist 


D 25 


L B 


logioLx 




ObsID 


Instr. 


Date 


Exposure 






(Mpc) 


(') 


(1O 1O L ) (lof 


;io(erg s -1 )) 


(10 20 cm- 2 ) 






(dd/mm/yy) 


(ks) 


High-Lx gal 


axies 




















IC 4296 


E Radio gal 


50.8^ 


3.8 


12.1 


41.54 


4.1 


3394 


S 


10/09/01 


25 


INLrO OU/ 


oA(r]0 


82. 6 3 


3.2 


14.9 


43.03 


5.4 


2882 


T 
1 


nc /m /no 
08/01/02 


43 


NGC 741 


E0 


75.8 3 


2.9 


14.2 


42.32 


4.4 


2223 


S 


28/01/01 


30 


NGC 1132 


E 


98. 2 4 


2.1 


8.3 


42.76 


5.2 


801 


S 


10/12/99 


13 


NGC 1399 


cD;El pec 


18. 5 1 


6.9 


4.2 


42.09 


1.3 


319 


s 


18/01/00 


56 
















4174 


1 


28/05/03 


45 


NGC 1600 


E3 


57.4 3 


3.2 


9.8 


42.19 


4.8 


4283 


s 


18/09/02 


22 


NGC 4472 


E2/S0(2) Sy2 


15. 1 1 


9.7 


7.5 


41.52 


1.7 


321 


s 


12/06/00 


32 
















322 


1 


19/03/00 


10 


NGC 5846 


E0-1;LINER HII 


21. 1 1 


3.8 


3.1 


41.57 


4.3 


788 


s 


24/05/00 


23 


NGC 7619 


E 


49.2 1 


2.6 


6.9 


42.06 


5.0 


3955 


s 


24/09/03 


31 


NGC 7626 


E pec 


51. 2 3 


2.7 


6.9 


41.51 


5.0 


2074 


I 


20/08/01 


26 


Moderatc-Lx galaxies 


NGC 720 


E5 


25. 7 1 


4.6 


3.1 


41.33 


1.5 


492 


s 


12/10/00 


29 


NGC 1332 


S(s)0 


21.3 1 


4.1 


2.3 


40.95 


2.2 


4372 


s 


19/09/02 


45 


NGC 1387 


SAB(s)0 


18.9 1 


3.3 


1.2 


40.78 


1.3 


4168 


1 


20/05/03 


45 


NGC 1407 


E0 


26.8 1 


5.3 


6.5 


41.32 


5.4 


791 


s 


16/08/00 


40 


NGC 1549 


E0-1 


18.3 1 


4.7 


2.9 


40.66 


1.5 


2077 


s 


08/11/00 


22 


NGC 1553 


SA(rl)0 LINER 


17.2 1 


5.3 


3.7 


40.64 


1.5 


783 


s 


02/01/00 


11 


NGC 1700 


E4 


41.1 1 f 


3.0 


1.6 


41.20 


4.8 


2069 


s 


03/11/00 


27 


NGC 3607 


SA(s)0 


21.2 1 


4.5 


3.4 


40.94 


1.5 


2073 


1 


12/06/01 


38 


NGC 3923 


E4-5 


21.3 1 


6.4 


4.9 


41.03 


1.5 


1563 


s 


14/06/01 


8.8 


NGC 4365 


E3 


19.0 1 


5.8 


3.7 


40.83 


6.2 


2015 


s 


02/06/01 


40 


NGC 4552 


E;LINER HII 


14.3 1 


5.0 


2.0 


40.65 


2.6 


2072 


s 


22/04/01 


54 


NGC 5018 


E3 


42. 6 3 


3.6 


7.1 


40.96 


7.0 


2070 


s 


14/04/01 


28 


Low-Lx galaxies 


NGC 3115 


SO 


9.0 1 


7.3 


1.4 


39.86 


4.3 


2040 


s 


14/06/01 


36 


NGC 3585 


E7/S0 


18.6 1 


6.1 


3.3 


40.35 


5.6 


2078 


s 


03/06/01 


35 


NGC 3608 


E2 LINER 


21.3 1 


3.2 


1.7 


40.39 


1.5 


2073 


1 


12/06/01 


38 


NGC 4494 


El-2 LINER 


15.8 1 


4.5 


2.3 


40.36 


1.5 


2079 


s 


05/08/01 


15 


NGC 4621 


E5 


17.0 1 


5.0 


2.5 


40.14 


2.2 


2068 


s 


01/08/01 


25 


NGC 5845 


E 


24.0 1 


0.9 


0.45 


40.17 


4.3 


4009 


s 


03/01/03 


30 



daxiq 
SBF: 



cam p 



Icnscn et al 



istanccs were 



(2003) 



od fr pm. 



D n -£r: Faber et al 



SBF: Tonry et al 



Note. — Listed above are all of the ; 
new Cepheid zero-point (see text), 

Lb was determined from the face-an_ reddening-corrected B-band magnitude given by LED A. 
consistently in the present work(§ p. 2j) . excluding obvious emission from any low-luminosity AGN, 
■fiducial 300 kpc radius. The galaxy type was taken from NED. Nu is the nominal Galactic column-density along the line-of-sight. We show the 



(2001), corrected for the the 
rcdshift distance (LEDA)\ t — uncertain. 

10.0 keV band self- 



fll989| ), 

Lx was computed in the 0.1- 

and extrapolating the surface brightness to a 



Chandra observation identifier (ObsID), the ACIS instrument (I or S) and the net cxposurc-timc, having excluded periods of flaring. 



account for it in our spectral fitting. 

3.1. Background estimation 

One of the key challenges in spectral-fitting diffuse X- 
ray emission is ensuring proper background subtraction. 
For Chandra a set of blank-field event files h ave been 
made available as part of the standard Caldl distribu- 
tion, from which background spectra can be accumulated 
corresponding to similar regions of the detector. For each 
observation, we prepared from these a suitably projected 
background events file. We were able to extract from this 
file "template" background spectra for each region of the 
detector in which our "source" spectra were obtained. 
However, these background spectra are unlikely to rep- 
resent perfectly the background in any one observation. 
There are known to be significant long-term secular vari- 
ations in the non X-ray components of the background, 
substantial field-to-field variation in the cosmic compo- 
nent, and there may be some residual mild flaring, ft is 
also worth noting that the hard (power law) component 
of the cosmic X-ray background arises from undetected 
background AGN, so its absolute normalization is also a 



strong function of the point source detection complete- 
ness; in turn this is a function of the sur face brightness of 
the ga laxy and the total exposure time ( Kim fc Fabbiano 
2004a|) . 

Several authors have adopted the practice of renor- 
malizing the background template to ensure good agree- 
ment with the instrumental background at high energies 
(> 10 keV). Such a procedure, however, also renormal- 
izes the (uncorrelated) cosmic X-ray background and in- 
strumental line features, which can lead to serious over 
or under-subtraction. Given these reservations we chose 
to use an alternative background estimation procedure. 
Our method involved modelli ng the backgr o und, some- 



(|2004|) . For 



what akin to the approach of Buote et al. 
each observation, we extracted a spectrum from a "source 
free" region of the ACIS field of view. We chose a ^2' 
region centred on the SI chip if the galaxy was cen- 
tred on S3, or on the S2 chip where the galaxy was 
centred on AC1S-1. If the SI or S2 chips were turned 
off, we chose a small, < 1' region on the S3 or ACIS- 
I chips, as appropriate, positioned to be in a region 
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of as low surface-brightness as possible. We excluded 
data from the vicinity of any point-sources found by the 
source detection algorithm. Additionally we extracted a 
"source+background" spectrum from the CCD on which 
the source was centred, in an annulus centred at the 
galaxy centroid and with an inner and outer radii typi- 
cally ~2.5' and 3.3'. We adopted two spectra since we 
found that this procedure enabled us most cleanly to 
constrain the background. Some of the brightest galax- 
ies are so extended in the X-ray that even in our "source 
free" region there is a small contribution from hot gas 
at large radii. We found that using two spectra with 
different hot gas contributions allowed this to be read- 
ily disentangled from the actual background components. 
In order to constrain the model, we fitted both spectra 
simult aneously, without background subtraction, using 
Xspet. Our model consisted of a single APEC plasma 
(to take account of the diffuse emission from the galaxy; 
the "source"), plus background components. These com- 
prised a power law with T = 1.41 (to account for the hard 
X-ray background), two APEC models with solar abun- 
dances and kT= 0.2 and 0.07 keV (to account for the soft 
X-ray background) and, to model the instrumental con- 
tribution, a broken power law model and two Gaussian 
lines with energies 1.7 and 2.1 keV and negligible intrin- 
sic widths. We have found that this model can be used to 
parameterize adequately the template background spec- 
tra. 

In order to disentangle the source and background 
components, given the general lack of photons in these 
spectra, we tied the abundances and temperatures of the 
"source" APEC components between both extraction re- 
gions, but allowed the normalizations to be free. In the 
fainter galaxies the normalization of the source compo- 
nent in the "source free" region tended, as expected, to 
zero. We also assumed that the background model nor- 
malization scales exactly with the extraction area. In 
general, we found that this was able to fit both spectra 
very well. In our subsequent spectral analysis, we did not 
background-subtract the data using the standard tem- 
plates, but took into account the background by using 
an appropriately scaled version of this model. 

Even if the background spectrum varies substantially 
over the field-of-view, our background modelling is most 
correct at largest distances from the source centroid, 
where the results are most sensitive to the background. 
In fact, we found that the standard background tem- 
plates fared much worse than these modelled background 
estimates when the data were from regions of low surface 
brightness. We discuss this issue further, and how the 
choice of background can affect our results in 



5.2 



3.2. Lx estimation 

In order to provide a self-consistent analysis of the 
galaxies in th is sample , we obtained estimates of Lx 
based on the Chandrc data. First, we estimated the 



flux within our chosen spectral extract ion regi ons from 
our the best-fitting spectral models (§ 4.3-4.4). Fluxes 
were computed separately for the gas and undetected 
point-sources in the energy-band 0.1-10.0 keV. Since the 
adopted aperture will not contain all of the diffuse flux 
from the galaxy, we extrapolated the emission out to a 
projected radius of 300 kpc, i.e. the Virial radius for a 
1.6 x 10 12 M Q galaxy. We assumed spherical symmetry 



and parameterized the surface brightness with a single 
or double /3-model. The /3-model parameters were de- 
termined from fits to the radial surface brightness in the 
0.3-2.0 keV band, using dedicated software which can 
fold in the instrumental point-spread function, which we 
computed at 1 keV. Data from the vicinity of any de- 
tected point-sources were excluded from the fit, and we 
assumed that the hot gas and undetected sources had the 
same radial brightness distribution. The results of the 
surface brightness fits to each galaxy will be discussed in 
detail in a subsequent paper. This procedure typically 
corrects the flux upwards by a factor ~1.1— 4, depending 
on the shape of the surface brightness profile. Since it is 
by no means certain that this extrapolation is valid out 
to ^300 kpc, we expect this to introduce some uncer- 
tainty into the estimated Lx- However, for our present 
purposes we believe this approach is sufficiently accurate. 

To compute the total Lx of the galaxy, we also included 
the flux of all detected point-sources within the B-band 
twenty-fifth magnitude (D25) isophote, all of which were 
assumed to be associated with the galaxy for these pur- 
poses. We fitted the composite spectrum of all these 
sources with our canonical (kT=7.3 keV bremsstrahlung) 
model. In most cases this gave a good fit to the data, 
although in a few instances a single power law or power 
law plus disk blackbody components were used instead 
to obtain a good fit. In the event a significant low- 
luminosity AGN appears to be present in the galaxy (i.e. 
NGC 1553, IC4296), we omitted the flux from the AGN. 
The total luminosity of the (detected, plus undetected) 
point-sources wi thin D25 was typically in agreement with 
the estimate of Kim & Fabbiano ( 2004a| ), based on ex- 
trapolating the resolved X-ray luminosity functions in 
nearby early-type galaxies. We found a mean Lx (point 
sources)/Le^ 0.9 x 10 30 erg s _1 L0 _1 , in excellent agree- 
ment with these authors. The point-source populations 
will be discussed in detail in a subsequ ent paper 



Com paring with the fluxes given in O 'Sullivan et al 



(2001) we find broad agreement, although our estimates 
tend to be ~0.25 dex higher. We attribute this discrep- 
ancy to differences in spectral modelling and our extrap- 
olation procedure. 

4. SPECTRAL ANALYSIS 

4.1. Solar abundances standard 



Throughout this paper 
spheric abundances of Asplund et al 



we adopt the s olar photo- 
" Q2004P , which de- 



viate significantly fo r many of the key species from the 
previous standard of Grevesse & Sauva] ( 1998| ). The in- 
corporation of detailed 3D line transfer modelling (and, 
in some cases, treatment of non-LTE conditions) has 
tended to reconcile so-called "meteoritic" and photo- 
spheric abundances (for non- volatile species) , so that dis- 
crepancies remaining are typically at <0.1 dex, i.e. ap- 
proximately the same level as the statistical uncertain- 
ties. We therefore adopt these abundances as our stan- 
dard. This does, however, introduce some complications 
when comparing with results reported in the literature 
since most authors adopt eithe r the older abundances 
standard of Grevesse fc Sauva| (1998 ) or the outdated 
abundances oF Anders & Grevesse "(1989). 

For comparison with our work, Zp c , Zo/Zp e , Zn /Zfc, 
ZMg/ZFoj Zs j/ZF n , Zr/Zf^ and Z^/Z^y , referenced to the 



standard of Grevesse & Sauval (1998) should be scaled 
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shown along with the quality of fit. Statistical errors represent the 90% confidence region. Fig 
an estimate of the sensitivity of each measurement to possible sources of systematic error (see § 



for all galaxies in the sample, calibration uncertainties (§ 



3 given in square brackets are 
Since we did not assess them 



are not included in this error-budget for most of the sample. These 



slit 

(§ ksD 



certainly not be added in quadrature with the statistical errors. Where we were able to obtain an abundance gradient 
we estimated an e mission- weighted Zpv. ex trapolated over a large aperture (see text); those affected galaxies are marked 



(^). — results taken from 
be constrained, they were fixed at the Solar 



Humphrey et al 



(2004). corrected to our adopted abundance standard. Where parameters could not 
value, and listed as (...)■ 
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by 1.12, 1.32, 1.55, 1.00, 0.98, 1.38 and 0.93, respectively. 
Likewise, th ese abundances obt a ined w ith the outdated 
standard of Anders fe Grevessej ( 1989 ) should be scaled 
by 1.66, 1.12, 1.07, 0.68, 0.66, 0.71 and 0.63, respectively. 

4.2. Error-bar estimation 

It is w orth n oting that the error-bars generated sim- 
ply with Xspec tend to underestimate substantially the 
true confidence region. This seems to arise because the 



fitting algorithms employed often spuriously identify fit 
convergence (regardless of how low one sets the "criti- 
cal delta-fit statistic" parameter). This problem can be 
mitigated by re-starting the fit a number of times from 
the apparent minimum and ensuring that the fit statis- 
tic has truly minimized. Unfortunately the error-bar es- 
timation routines assume that a single iteration of the 
fitting algorithm is sufficient to characterize the shape of 
the x 2_ space, which is seldom true. Since the failure to 
converge results in an inflated \ 2 f° r anv given value of 
an interesting parameter, this produces error-bars which 
are systematically too small, sometimes dramatically so. 

We experimented with two strategies to overcome this 
problem, which we found to agree well in general. In 
our preferred tec hnique , we used a script to emulate the 
behaviour of the Xspec "error" command, but minimiz- 
ing x 2 a t each step with multiple iterations of the fit 
command (until the \ 2 converged within a tolerance of 
0.001). This has the advantage of identifying (small) 
local minima. The alternative met hod, extensively cm - 
ployed in our previous papers (e.g. Buote et al. 20031) 
involves performing Monte-Carlo simulations in which 
spectra are simulated from the best-fit model. The best- 
fit model is then fitted to each simulated dataset and the 
confidence region inferred from the distribution of the 
best-fitting parameter values. Provided the fit statistic 
is well-described by the x 2 distribution, both techniques 
are statistically equivalent. 

4.3. Spatially-resolved spectroscopy 

Where the data were of sufficient quality, we extracted 
spectra in a number of concentric annuli, centred on the 
nominal X-ray centroid. We determined the centroid it- 
eratively by placing a 0.5' radius aperture at the nominal 
galaxy position (obtained from NED) and computing the 
X-ray centroid within it. The aperture was moved to the 
newly-computed centroid, and the procedure repeated 
until the computed position converged. Typically the X- 
ray centroid agreed with that from NED. The widths of 
the annuli were chosen so as to contain approximately the 
same number of background-subtracted photons and en- 
sure there were sufficient photons in each to perform use- 
ful spectral-fitting. We restricted the fits to the energy- 
band 0.5-7.0 keV, to minimize instrumental background, 
which dominates at high energies, and to avoid calibra- 
tion uncertainties at lower energies (however, see § |5.4|) . 
The spectra were rebinned to ensure a signal-to-noise ra- 
tio of at least 3 and at minimum 20 photons per bin (to 
validate x 2 -fitting). 



We fitted the spectra using [Xspeq with a model com- 
prising a hot gas (vapec) component, plus an addi- 
tional 7.3 keV bremsstrahlung component to take ac- 
count of undetected point-source emission (this model 
give s a good fit to th e detected sources in nearby galax- 
ies: Irwin ct al. 2003). We used a slightly modified form 



of the existing Xspec] vapec implementation so that Zp e 
is determined directly, but for the remaining elements 
the fit parameters were the abundance ratios (in solar 
units) with respect to Fe. This model enabled errors 
on the abundance ratios to be determined directly from 
Xspec. In general, the data did not enable us to con- 
strain any abundance ratio gradients, and so we tied the 
abundance ratios between all annuli. The absorbi ng col- 



umn density (7Vh) was fixed at the Galactic value (Dickey 



fe Lqckman |199C ) ; the effect of varying Nr is discussed 
in § [5l| Where abundances could not be constrained, 
they were fixed at the Solar value. We took account of 
possible hot gas spectral projection effects by employing 
the projet model implemented in Xspec, where possi- 
ble. By fitting the data in several annuli, it is possible 
to measure and, crucially, constrain, any temperature 
gradients which contribute to the so-called Fe bias. We 
discuss the temperature profiles of each galax y in detail 
in a subsequent paper, and briefly in § [4.3.1 . In order 



to improve the abundance constraints, which tend to be 
somewhat poorer than the temperature constraints, it 
was sometimes necessary to tie the abundances between 
adjacent annuli (see below). In the interests of physi- 
cally meaningful results, we constrained all abundances 
and abundance ratios to the range 0.0-5.0 times the so- 
lar values. Abundance ratios were kept fixed at the solar 
value unless we could obtain interesting constraints dur- 
ing the fitting. The spectra for the central annuli of all 
the galaxies are shown in Fig [l| — along with the best- 
fitting model. Table|| lists globally-averaged abundances 
derived from each fit. 

Since there is some evidence of limi ted multi-phase gas 
in some gia n t ellipticals and groups ( Buote et al. 2003a ; 
Buote 2002 ; Xue et al. 2004 ), we experimented with the 
addition of an additional hot gas component to the inner 
few annuli of each galaxy. In a few cases (see below) this 
improved the fit significantly. It was generally not pos- 
sible to determine the abundances of this additional hot 
gas component separately, and so they were tied to those 
of the other gas component in the same annulus. This 
two component model is a simple parameterization of 
multiple temperature gas components in the extraction 
aperture, for example if there is a strong temperature 
gradient over t he extraction regio n. It has been shown 
previously (e.g. Buote et al. [2003b| ) that abundance con- 
straints obtained with this parameterization are gener- 
ally in very good agreement with those derived from 
more complex models which allow the emission-measure 
to vary continously as a function of temperature. In all 
cases discussed below we included a component to ac- 
count for undetected point-source emission. 

To estimate an emission-weighted global abundance 
where there was some evidence of an abundance gradient, 
we first parameterised the radial surface brightness pro- 
file of each galaxy with a /3-modcl fit. We then assumed 
that the underlying abundance profile is well-described 
by a broken power law model. This model is a good fit 
to the profile seen in NGC 1399, and is similar to that in 
the centre of the group NGC 5044 (|Buote et al.| |2003b| ). 
Where there was insufficient data to constrain the ra- 
dius of the break, we simply fixed it to the midpoint of 
the innermost data-bin, inside which the abundance is 
flat (to prevent it inflating unphysically at small radii). 
Using this model we were able to estimate the emission- 
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Fig. 1. — Spectra of the high-Lx galaxies. Where spatially-resolved spectroscopy was possible, only the central bin is shown, for clarity. 
In addition we show the best-fit model folded through the instrumental response (solid line; red), and the contribution from undetected 
point sources (dash dot line; green) and the one or two hot-gas components (dashed lines; blue and magenta). The energies of interesting 
line features are included to guide the eye. 



weighted abundance extrapolated to a "global" aperture 
(for which we used 30'; the results are relatively insen- 
sitive to this choice). For those galaxies with measured 
abundance profiles, we show the data and the best-fit 
model in Fig || 

4.3.1. Comments on individual galaxies 

NGC 507. We extracted spectra in seven contiguous, 
concentric annuli, with outer radii 0.6, 1.0, 1.3, 1.9, 2.7, 
4.0 and 5.6' (14, 23, 31, 44, 63, 93 and 130 kpc), respec- 
tively. To improve the abundance constraints, we tied 
Z Fc between the 1 st and 2 nd annuli, the 3 rd and 4 th an- 
nuli and the 5 th and 6 th annuli. The temperatures were 



allow ed to vary in each annulus. The best-fitting temper 
ature profile rose trom ~1 keV in the centre to ~1.3 keV 
in the outermost radii. The abundance profile is essen- 
tially flat, except for the outermost bin. The m e asure d 
abundances generally agreed with Kraft et al. (2004), 



who used the Chandra data. We did not find any statis- 
tically significant improvement in x 2 if a second hot gas 
component was added (Ax 2 ~2). Nevertheless, to effect 



a comp arison with the XMM results of Kim fc Fabbiano 
(2004b), we experimented with the addition of such a 
component to the innermost bin. We found that this in- 
creased the central Zp e to 0.8 7^o Js, which is marginally 
inconsistent with Zp c rvj 2 found by these authors. The 
reason for this discrepancy is unclear since the temper- 
atures of these components (kT~0.8 keV and 1.4 keV) 
and the total flux of unresolved point-sources within D25 
(~ 4 x 10~ 13 erg s^'em" 2 ), which might system atically 



affect Zfp . were v ery similar to those found by Kim 
Fabbiano (|2004b|) . Nonetheless, our abundance ratios 



were broadly consistent with those reported for XMM. 

NGC 741. We extracted spectra in three contiguous, 
concentric annuli, with outer radii 0.6, 2.0 and 3.8' (12, 
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43 and 82 kpc), respectively. In order to obtain inter- 
esting constraints, we tied Zp between the inner two 
annuli. We found a statistically significant (Ax 2 = 10) 
improvement in the fit if two, rather than one, hot gas 
components (kT=0.67 keV and 1.2 keV) were used in the 
innermost bin; there is no evidence for the cooler compo- 
nent in the outer bins. The best-fitting abundances are 



somewhat higher than Zf c =0.371o is found by Mulchaey 



et al. Q2003] ), who fitted Rosat data with a single MEKAL 
model (and no undetected source component). We at- 
tribute the discrepancy to the Fe bias. 

NGC 1132. We extracted spectra in four contiguous, 
concentric annuli, with outer radii 0.8, 1.6, 2.4 and 3.7' 
(21, 43, 66 and 100 kpc), respectively. In order to ob- 
tain interesting constraints, we tied the abundances in 
all four annuli. We found a ~flat temperature profile 
(kT~1.0 keV) when only one hot gas component was 
used. We found a slight improvement in the fit if two hot 



gas components (kT=0.8 keV and 1.6 keV) were used in 
the central bin. Our best-fitting abundance within the 4' 
radius was in good agreement with previous ASCA de - 



199^ ) 



terminations (Buote 2000b; Mulchaey fc Zabludofi] . 
Recent XMM measurement s show evidence of a sig mf- 
icant abundance gradient flGastaldcllo ct aL 2004| , F. 



Gastaldello et al. 2005, in prep.); within ~4', our abun- 
dance determination is in good agreement with XMM. 

NGC 1399. We extracted spectra in 10 contiguous, 
concentric annuli, with outer radii 0.2, 0.7, 1.5, 2.4, 3.0, 
3.6, 4.2, 5.4, 7.7 and 12' (1, 4, 8, 12, 16, 20, 22, 29, 41 
and 64 kpc), respectively. In order to improve abundance 
constraints, we tied together Zp c between annuli 2 and 3, 
between annuli 4-6 and between annuli 7-9. We found a 
significant improvement in the fit if two hot gas compo- 
nents were used in the inner 6 bins. For the cooler compo- 
nent kT rises from ~0.7 keV to 1.3 keV. The hotter com- 
ponent has kT^1.5 keV, but was less well-constrained. 
Although the best-fitting model was not formally accept- 
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Fig. 3. — Spectra of the low-Lx galaxies. 
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Fig. 4. — The Fe abundance profiles, for those systems in wl 
models used to infer a global emission-weighted abundance (§ 
show separately for NGC 1399 the profile obtained when two ho' 



ich an abundance gradient was measured. Also shown are the best-fitting 
1.3 1. To demonstrate the Fe bias and the importance of deprojection, we 
;as components were used in those annuli requiring them (solid lines) and 



when only a single hot gas component is used (dotted lines). We also show results for which deprojection was used (marked with squares) 
and for which deprojection was not (marked with triangles). 



able adding an additional hot gas component did not im- 
prove the fit further. Given the excellent S/N of the data 
and the fact that the fractional fit residuals are typically 
~ a few percent, it seems probable that the remaining 
errors are primarily systematic, for example calibration 
effects. We found evidence of an abundance gradient, 
with a sharp drop-off at ~7', as shown in Fig |[ Our 
abundance profile is in e xcellen t agre ement with that de- 
rived from XMM data ( Buotej 2003) . These results are 
also in agreement with previo us single- aperture AS CA 
abundance measurements (e.g. Buote & Fabian 1998). 

NGC 1407. We extracted spectra in 3 contiguous, con- 
centric annuli with outer radii of 0.3, 0.8 and 2.1' (2.7, 6.1 



and 16 kpc), respectively. In order to improve the con- 
straints, the abundances were tied between the two inner 
annuli. We only required a single hot gas component in 
each annulus. The temperature rose from ^0.65 keV in 
the centre to 1.0 keV in the outer annuli. We find some 
evidence of an abundance gradient in the data, although 
the error-bars are rather large. Given the large errors on 
our interpolated global abundance, our best-fitting val- 
ues are in general agreement with the ASCA value of 
Z Fc = 0.88^60 found by |Buote fc Fabia"^ fll998| ). 



NGC 1600. We extracted spectra in 3 contiguous, con- 
centric annuli with outer radii 0.8, 2.2 and 3.8' (13, 
36 and 62 kpc), respectively. In order to obtain inter- 
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esting constraints, we tied Zp c between the inner two 
annuli. We found a significant improvement in the fit 
(A% 2 =15) when two hot gas components (kT=0.86 keV 
and ^3 keV), rather than one, were used in the central 
bin. In the outer radii, only a sin gle hot g as co mponent 
(kT=1 .5 keV ) was needed. Using Chandrc data Sivakofi 
et al| ( 2004 ) reported Zp ~1.8 for a two-temperature fit 
within 1 effective radius, falling to ^0.5 at ~3' (using 
our solar standard), which is consistent with our results. 

NGC 4472- We extracted spectra in 9 contiguous, con- 
centric annuli with outer radii 0.3, 0.6, 1.1, 1.6, 2.2, 3.0, 
4.0, 5.2 and 7.7' (1.3, 2.8, 5.0, 7.1, 9.6, 13, 17, 23, 34 kpc), 
respectively. In order to improve the constraints, we tied 
together Zp e in annuli 2 and 3, and also in annuli 4-7, 
and annuli 8-9. Only a single hot gas component was re- 
quired by the data in any given radius. The temperature 
profile rose from kT=0.67 keV in the centre to 1.3 keV 
in the outer annuli. There is some evidence of slight 
abundance gradient, although the error-bars are rather 
large. Our measured Zp c in each bin is somewhat larger 
than the values ~0.5-l.l (converting to o ur solar stan- 
dard) found by Finoguenov fe Jones ( 2000 ). This is most 
likely a consequence of the Fe bias; there is a clear tem- 
perature gradient even within the smallest annuli those 
authors a ttempt to use and furth ermore our results agree 
well with |Buote fc Fabiaij ( |l998|) . 

NGC 4552. We extracted spectra in 3 contiguous, con- 
centric annuli with outer radii 0.3, 1.2 and 2.4' (1.2, 
5.2 and 10 kpc), respectively. We tied abundances be- 
tween all annuli. Only one hot gas component was re- 
quired in any annulus to fit the data. There is some 
evidence of a central temperature peak, sim ilar to that 
found in NGC 1332 (|Humphrey et al.| |2004| ): we found 
kT=0.58 keV in the centre, falling to ~0.33 keV by the 
next annulus. Using ASCA Finoguenov & Jones ( |200C| ) 



found Zp c < 0.12, when fitting the spectrum with a sin- 
gle hot gas component. We attribute this discrepancy 
to the Fe bias. The abundance ratios reported by those 
authors also disagree with our measured values, which is 
presumably a related effect. 

NGC 5846. We extracted spectra in 6 contiguous, con- 
centric annuli with outer radii 0.4, 0.7, 1.7, 2.3, 3.9' (2, 
4, 7, 10, 14 and 24 kpc), respectively. We found that fit- 
ting a single hot gas component (plus undetected source 
component), with kT rising from ~0.6 keV in the cen- 
tre to ~1 keV in the outer radii, gave a rather poor fit 
(x 2 /dof=630/450). Adding an additional hot gas compo- 
nent to the inner annuli did not improve the fit. However, 
the fit was considerably improved if such a component 
was added to the outermost annulus (A% 2 =60). By in- 
spection, the image shows considerable large scale struc- 
ture, which co mplements the disturbed mor phology on 
smaller scales ( Trinchieri fc Goudfrooij 2002 ). To inves- 
tigate this further, we constructed a hardness map (using 
energy-bands 0.1-0.8 keV and 0.8-7.0 keV, separately 
smoothing each image by convolution with a Gaussian 
of width 5"), which revealed that the gas temperature 
within each annulus is not uniform. We confirmed this 
explicitly for the outermost annulus by restricting the ex- 
traction region to a narrow sector chosen to contain the 
softest photons; the best-fitting temperature of the gas 
in the outer annulus was then significantly lower (falling 
to 0.8 keV). A detailed analysis of the two-dimensional 
temperature structure of this galaxy is beyond the scope 



of this present work. However it is sufficient to add ad- 
ditional hot gas components in each annulus for which 
the fit was significantly improved to obtain a reasonable 
estimate of Zp - The addition of a third hot gas compo- 
nent to each annulus did not improve the fit significantly. 
Although it is not clear that spherically-symmetric de- 
projection is appropriate in such a system, the fit was 
somewhat better when projct was used than when it was 
omitted. The final \ 2 was poor; this may reflect system- 
atic errors in the response, since the source is bright, or it 
may reflect the imperfect spectral dcprojcction resulting 
from the spherical approximation, or it may be a con- 
sequence of the complicated temperature structure. We 
did not find any clear evidence of an abundance gradient 
so the abundances of all gas components were tied. 

NGC 1619. We extracted spectra in 7 contiguous, con- 
centric annuli with outer radii 0.2, 0.6, 1.1, 1.8, 2.5, 3.4 
and 4.5' (2.5, 8.0, 16, 25, 35, 48 and 63 kpc), respectively. 
We did not find strong evidence of a Zp G abundance gra- 
dient and so we tied abundances between all annuli. Only 
one hot gas component was required in each annulus to fit 
the data. The temperature was seen to rise from ^0.76 
keV in the centre to ~0.95 keV in the outer annulus. 
Our measur ed Zp c is in exc e llent agreement with that 
reported by |Buotc fc Fabiarj ( |1998[ ). 



4.4. Single aperture spectroscopy 

In many cases there were insufficient photons to at- 
tempt spatially-resolved spectroscopy. We therefore ex- 
tracted the spectra from a single aperture, the size of 
which was chosen crudely to maximize the S/N, and 
was typically ~2-3'. We fitted our canonical hot gas 
plus bremsstrahlung model to each spectrum. We note 
that, for low-Lx systems, this model is tantamo unt to 
the two-tempera ture hot gas modelling used by Buotc 
Fabian ( 1998 ), who found that one component was 
typically hot (kT> 5 keV) in the most X-ray faint sys- 
tems. In some cases a statistically significant improve- 
ment in the fit statistic (< 1% chance the improvement 
is spurious, on the basis of an f-test) was found if an 
additional hot gas component was included. We state 
below, for each galaxy, whether one or two hot gas com- 
ponents were required. All fits included an undetected 
point-source component. The spectra for all the galaxies 
are shown in Fig ]l]-[| along with the best-fitting model. 

4.4.1. Comments on individual galaxies 

IC4296. We extracted the spectrum from a 1.5' 
(22 kpc) aperture, omitting data from the vicinity of the 
central LLAGN. We found that the addition of an extra 
hot gas component significantly (A% 2 = 11) improved the 
fit. The gas components had kT=0.7 keV and 1.3 keV. 
Our results a re in excellent agreem ent with Zp e derived 
from ASCA Quote fc Fabian||l998| ). 

NGC 1387. We extracted the spectrum from a 2' 
(11 kpc) aperture. Only one hot gas component 
kT=0.65 keV) was needed. Using Rosat lones et al 



1997|) measured ZF e =0.25^Q22, in agreement with our 
results, even though these authors omitted a component 
to account for point-source emission, which contributes 
~35% of the flux in our extraction aperture. 

NGC 1549. We extracted the spectrum from a 3' 
(16 kpc) aperture. There was very little diffuse emission, 
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but some evidence of metal enri chment in the ~0 - 4 keV 
gas. Using the Rosat PSPC, [Davis fc White] fll996| ) 



reported Zp e < 0.18, omitting a component to account 
for undetected point sources. Point sources contribute 
~50% of the flux in our extraction region, and so we at- 
tribute the discrepancy to incorrect modelling of these 
sources, which could not be resolved with Rosat. 

NGC 1553. We extracted the spectrum from a 2' 
(10 kpc) aperture, e xcluding d ata in the vicinity of the 
central LLAGN. The Chandra image reveals remarkable 
structure (a striking S-shaped pattern, wh ich may be re- 



lated to interaction with the central AGN: Blanton et al 



Eooip 



Blanton ct al.| fl200l|) rep orted Z Fc = OTTTf^ for 



the hot gas, using the Chandrc data, which agrees very 
well with our fit with a single ^0.4 keV gas component. 
We do not find any evidence of a hard component, such 
as reported by these authors, in addition to the hot gas 
and point source components. However, if we include a 
small amount of flaring-contaminated data, such a com- 
ponent is required. Given the low abundance and the 
complexity of the ISM, we experimented with the addi- 
tion of a second (kT~0.86 keV) gas component. It did 
not improve the fit appreciably (A% 2 ~ 2), but Zp c was 
significantly increased to 0.24+°;^, and 2^ e /Zp e reduced 
to O.67I0 57. The very low abundance we obtain for this 
galaxy may, therefore, be mitigated by the Fe bias. How- 
ever the data do not require such an additional hot gas 
component. 

NGC 1700. We extracted the spectrum from a 1.5' 
(18 kpc) aperture. Only a single ~0.4 ke V gas co mpo- 
nent was required to fi t the data. Using Chandra data 
' tatler fc McNamara| ( |2002p reported Z Fc =0.83±0.50, 
omitting a component to account for undetected point- 
sources (which contribute~30% of the flux in our extrac- 
tion region). Nonetheless, this value is in agreement with 
our results. 

NGC 3115. We extracted the spectrum from a 2' 
(5 kpc) aperture. The diffuse component was dominated 
by the undetected point-source component. However, 
there is still some evidence of a single ~0.4 keV gas com- 
ponent. The abundances were very poorly constrained. 

NGC 3585 and NGC 4494. We extracted spectra for 
both of these sources from 2' (11 and 9 kpc, respectively) 
apertures. Both galax i es we re studied with XMM by 



i 



ICTSullivan fc Ponman| (|2004j), who reported Z Fc < 0.1, 
with tight abundance constraints. Such a low abun- 
dance is excluded by our data for NGC 3585, although 
it is consistent within the very lar ge error-b ars we found 
for NGC 4494. Inspection of the Chandra spectrum of 
NGC 3585 clearly reveals the presence of the Fe hump, 
indicating substantial Fe enrichment (Fig H). In con- 
trast the hot gas component in NGC 4494 is clearly over- 
whelmed by the point-source contribution. The lower 
spatial resolution of XMM makes s ubstantial ly more 
point-source contamination than for 
spectra inevitable, 
with our results, since the metal abundance of the hot 
gas will then be highly sensitive to the ability to model 
the spectrum of the undetected sources. Although on 
average the composite spectra of detected sources in 
\Jhandr\ fields can be well-approxi mated by a simpl e 
bremsstrahlung or power law model (Irwin et al. 2003), 
there is no reason to believe that this is the exact spectral 



Chandra in both 
This may explain the discrepancy 



model for the undetected sources seen in any given galaxy 
(see § 5J)). This is problematical where the gas param- 
eters depend sensitively on its shape. Since Z Fe < 0.1 is 
difficult to understand in terms of the current picture of 
metal enrichment, but is more consistent with what we 
might expect if undetected sources are not properly ac- 
counted for, it seems likely that our best-fitting values 
are more representative of the true abundances in these 
systems. 

NGC 3601. A spectrum was extracted from a 2' 
(12 kpc) aperture. There is significant diffuse gas, clearly 
showing the characteristic Fe L-shell "hump" . The data 
only required one kT=0.46 keV gas component, and Z Fc 
is moderately well-constrained. Matsushita et al. (200C) 
reported a similar Z Fe using ASCA data. 

NGC 3608. A spectrum was extracted from a 2' 
(12 kpc) aperture. Only one ~0.35 keV gas component 
was required to fit the data. The hot gas contributes 
^50% of the flux in the extraction aperture, but the 
abundances were poorly constrained. 

NGC 3923. A spectrum was extracted from a 1.5' 
(9 kpc) aperture. We restricted spectral-fitting to the 
0.6-4.0 keV band, since there was some evidence of fea- 
tures outside this range which may be an artefact of 
the prolonged mild flaring in this observation, which it 
was impossible entirely to excise. Only one ~0.38 keV 
gas component was needed, in addition to undetected 
sources. The best-fitting Z Fo was poorly-constrained, 
but in good agreement with previous ASCA measure- 
ments ( Buote fc Fabian 1998 ). Intriguingly, the Zsi/Z Fc 
ratio was significantly higher than in the other systems, 
since producing the prominent Si line evident in the spec- 
trum (Fig ^ in such cool plasma requires an extremely 
high metallicity. This is not a very high significance ef- 
fect, since an acceptable fit (x 2 /dof=38/36) can still be 
obtained if we constrain Zgi/Z Fc = 1. An alternative ex- 
planation might be the presence of multiple temperature 
components in the aperture, in which case our inferred 
temperature may not be representative. There was some 
indication of an improved fit if an additional ^0.85 keV 
gas component was included, but the improvement was 
not significant (A% 2 =3). Fitting two temperatures did 
result in a smaller Si/Fe ratio (2.2^ 2 ' 2 ), which is large 
but marginally consistent with solar. 

NGC 4365. We extracted a spectrum from a 2' (11 kpc) 
aperture. We found two hot gas components (0.36 and 
0.93 keV) were required by the data, and the result 



ing Zp was poorly-constrained. Bivakoff et al. ( 2003| ) 



reported similarly p oorly-cons trained abundances for a 
similar region, using Chandr\ data. These authors also 
report some evidence of an abundance gradient. How- 
ever, given the challenges of background subtraction in 
such a low surface-brightness regime, and the few avail- 
able photons, we did not find compelling evidence. 

NGC '4621. We extracted a spectrum from a 2' (10 kpc) 
aperture. There is very little diffuse emission, and the 
undetected point source component dominates the spec- 
trum, however a ~0.3 keV gas component was needed. 
The abundances could not be constrained well. 

NGC 5018. We extracted a spectrum from a 2' (25 kpc) 
aperture. Only one ~0.5 keV gas component was re- 
quired to fit the data. The abundances were poorly con- 
strained. 
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NGC 5845. We extracted a spectrum from a 2' (14 kpc) 
aperture. There was very little diffuse emission, but some 
evidence of Fe-enriched ~0.83 keV gas. 

NGC 7626. We extracted a spectrum from a 2' (29 kpc) 
aperture. Only a single hot gas component, with 
kT=0.74 keV, was required. Our best-fitting Zp c were 
i n good agreement wit h those determined using ASCA 
flBuotc k Fabian||l998D . 



5. SYSTEMATIC ERRORS 

In this section, we address the extent to which sys- 
tematic uncertainties may impact upon our results. An 
estimate of the uncertainty due to these effects is given 
for each galaxy in Table 0. These numbers reflect the 
sensitivity of the best-fitting parameter values to each 
source of potential error, and we stress that they should 
certainly not be added in quadrature with the statistical 
errors. We discuss a number of effects in detail below. 
Those readers uninterested in the technical details of the 
analysis may like to proceed directly to § ^|. A break- 
down of the systematic error-budget for three represen- 
tative galaxies is shown in Table [| These are chosen 
approximately to span the luminosity, temperature and 
metallicity range of those galaxies in our sample for which 
interesting constraints were found. 

5.1. Calibration issues 



Uncertainties in the absolute calibration of Chandra 
may impact on our ability to determine reliable abun- 
dances for some or all of the sample. For example, it 
is not possible in standard processing to correct data 
from t he ACIS-S3 chip for "c harge transfer inefficiency" 
(CTI). Townsley et al. ( 2002[ ) proposed an algorithm to 
compensate for this effect, which we previously found to 
have a noticeable effect on Zp e of NGC 1332 (AZp e — 0.3) 
( Humphrey et al. 2004 ). Unfortunately it was not clear 
to what extent the change arose due to the effects of CTI 
or the older calibration of that technique. 

As an alte rnative m eans of assessing the calibration un- 
certainty of Chandrc , we instead reduced XMM archival 
observations of sever al galaxie s. Since XMM is calibrated 
independently from Chandrc\ , any discrepancies will be 
indicative of the magnitude of possible calibration er- 
rors. In the interests of efficiency, we only considered 
three galaxies for this analysis, NGC 1399, NGC 4552 
and NGC 3607, which reasonably span the temperature, 
metallicity and luminosity range in our sample . A fully 
self-consistent analysis combining XMM and Chandra 
data is beyond the scope of this present work, and will 
be addressed in a future paper. The XMM data were 
processed following a standard procedure outlined in F 
Gastaldello et al. (2005, in prep.) and spectra w ere ex- 



trac ted from similar apertures to those used in the Chan 
dra abundance analysis, where possible. 

For NGC 1399, we extracted spectra from 5 concentric 
annuli, with outer radii 1.5, 3.3, 5.5, 8.6 and 12 '. We fit- 
ted the same best-fitting model as used for the Chandra 



analysis. We obtained b est-fitting abundances in very 
good agreement with the Chandrc data, with the excep- 
tion of Zs/Zpo, which we found to be significantly lower 
(~0.65). 

For NGC 4552, we used a single 1.25' aperture, to 
which we fitted two hot gas com ponents (s ince a temper- 
ature gradient was evident in the Chandra data), plus an 



undetected source model. The temperatures of the com- 
ponents we re in goo d agreement with the range of kT 
resolved by Chandra . We found a significant reduction 
in the Zp e (AZp c — 0.3), comparable to that seen when 
CTI-correction was applied to the NGC 1332 data. We 
also found significant increases in Zo/Zpe, ZN e /ZFe and 
ZMg/ZFo, although these may be tied in to the reduction 
in Z Fe - 

For NGC 3607, we extracted data from a single 1.7' 
aperture (avoiding a chi p-gap) . W e found good agree- 
ment with our measured Chandra abundances, although 
the XMM constraints were somewhat poorer. The XMM 
gas temperature was slightly higher (kT~0.55 keV), al- 
though we foun d that th e temperature of this system 
determined with Chandrc was sensitive to the treatment 



of the background (more so than Zp c ), which probably 
explains this discrepancy. 

5.2. Background modelling 

One of the chief sources of systematic uncertainty in fit- 
ting extended, low surface-brightness objects is the treat- 
ment of the background. In our default analysis, we mod- 
ell ed the background, in a procedure akin to that adopted 
by Buote et al. ( 2004 ). We consider this to be the most 
robust approach, since it implicitly takes into account 
long-term variations in the non X-ray background (which 
are typically ~10%, in the absence of flaring 3 ), and field- 
to-ficld variations in the X-ray background. Further- 
more, it provides a natural means of subtracting off any 
extended, unrelated emission, such as cluster emission 
for galaxies in a cluster environment. 

To assess the sensitivity of our results to our back- 
ground treatment, we experimented with altering the 
background normalization by ±15%, comparable to the 
variation between background fields. We found little im- 
pact on our results for the highest surface-brightness re- 
gions, as might be expected. In contrast, however, the 
results for the faintest systems were rather sensitive to 
the background level. 



W e also experimented with using the standard Chan 
drc background templates. We found that they typically 



were not sufficiently accurate to be of use in very low 
surface-brightness regimes. For example, for the outer 
two annuli of NGC 1399, we found that using the tem- 
plates led to a substantial over-subtraction of the back- 
ground, giving a poorer fit (Ax 2 =128) and much higher 
abundances in these annuli (AZf o ^0.3-0.5 ), in disagree - 
ment with the XMM abundance profile (Buote 2002) ) . 
We reiterate that the common practice of renormalizing 
the background templates to match the high-energy (i.e. 
non-X ray background) is potentially dangerous since it 
will simultaneously renormalize the (uncorrelated) cos- 
mic X-ray background and instrumental lines. 

5.3. The plasma codes 

In order to estimate the extent to which the Fe L-shell 
modelling may impact upon our results, we experimented 
with replacing the (default) APEC plasma code with a 
MEKAL mo del. As has previously been seen (e.g. Buote 



(st alj|2003b| ), we found that this change led to slightly 
lower best-fitting Zp c where only a single temperature 
hot gas component was required in any annulus, and 
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TABLE 3 

Abundance measurement error budget 



Par. value AStat. Acalib Abkd Acode Aprojection Abandw. Asources A./Vjj AFe bias 



NGC 1399 


ZFe 


1.18 


±0.08 


-0.1 


±0.2 


±0.27 


±0.01 


+0.6 

-0.4 


±0.3 


-0.3 


-0.3 


Zo/Zp e 


0.41 


±0.06 


+0.02 


±0.02 


-0.11 


±0.01 


±0.13 


±0.03 


±0.03 


±0.03 


Znc/Zfo 


0.64 


±0.25 


-0.17 


±0.10 


-0.58 


-0.04 


+0.04 
-0.21 


+0.05 
-0.14 


-0.12 


-0.24 


ZMg/^Fc 


0.76 


±0.08 


+0.08 


+0.08 
-0.04 


-0.22 


+0.02 


-0.17 


±0.02 


-0.17 


±0.10 


Zsi/Zp c 


0.84 


±0.06 


+0.04 


±0.05 


-0.17 


±0.02 


-0.07 


±0.02 


-0.08 


+0.11 


Zs/Zf 


1.14 


±0.12 


-0.49 


-0.07 


-0.38 


±0.01 


±0.1 


±0.01 


-0.1 


+0.2 


Zni/Zfo 


2.61 


±0.32 


-0.86 


+0.4 
-0.1 


-1.0 


±0.1 


-0.2 


±0.2 


-0.3 


±0.7 


NGC 4552 


Zfo 


0.71 


+0.26 
-0.16 


-0.30 


+0.02 
-0.07 


-0.04 


±0.02 


+0.35 
-0.16 


-0.07 


±0.05 


±0.20 


Zo/ZFe 


0.17 


+0.11 
-0.06 


+0.11 


±0.03 


±0.07 


±0.01 


±0.21 


±0.02 


±0.01 


-0.03 


Znc/Zfo 


0.62 


±0.22 


+0.99 


±0.02 


±0.10 


±0.01 


+0.15 
-0.06 


±0.06 


±0.04 


±0.04 


ZMg/ZFo 


0.76 


±0.10 


+0.30 


+0.03 
-0.05 


±0.34 


±0.01 


+0.13 
-0.04 


±0.05 


±0.06 


±0.06 


Zsi/ZF C 


1.1 


±0.4 


+0.05 
-0.02 


±0.01 


±0.2 


±0.01 


+0.3 
-0.1 


±0.07 


±0.1 


±0.02 


ZNi/ZFc 


2.7 


±2.2 


-0.95 


±0.3 


-1.4 


±0.2 


±1.3 


-0.04 


±0.01 


-0.7 


NGC 3607 


Zfo 


0.32 


+U.60 
-0.14 


+0.08 


+ 11.1)4 
-0.11 


-0.06 




+ 1.2 

-0.11 


-0.06 


±0.010 


±0.06 


ZMg/ZFo 


0.63 


±0.63 


-0.33 


-0.29 


±0.28 




+0.04 
-0.26 


±0.09 


±0.02 


±0.06 



Note. — Emission-weighted abundance measurements for three representative galaxies chosen from the sample. In 
addition to the 90% statistical (Stat.) errors, we present an estimate of the possible magnitude of uncertainties on the 
abundances due to calibration (calib), background (bkd), plasma code (code), projection effects (projection), bandwidth 
(bandw.) and undetected source spectrum (sources). We also include an assessment of the impact of the "Fe bias" (AFe 
bias). For NGC 1399, this is the effect of fitting a only one gas component to each annulus (see text). For NGC 4552 and 
NGC 3607 it is the effect of adding an extra hot gas component to the central bin. It is difficult to assess the impact of 
multiple sources of systematic error simultaneously affecting our results (and it certainly is not correct simply to combine 
the errors in quadrature), so it is probably most correct to assume that the source of the largest systematic error dominates 
the systematic uncertainties. The magnitudes of the systematic errors listed are the changes in best-fit values, which can 
be more sensitive to such errors than the confidence regions. 



a slightly higher Zp e where two temperatures were re- 
quired. We found a complementary effect upon the abun- 
dance ratios, which may reflect in part the changing Zf c - 
These changes do not affect significantly the interpreta- 
tion of our results. Furthermore, we typically found that 
the MEKAL model gave slightly poorer fits to the data, 
justifying our preference for APEC. 

5.4. Bandwidth 

Presumably in part due to slight systematic errors 
in the response matrices, it is well-established that the 
choic e of bandwidt h affects the best-fitting abundances 
(e.g. Buote 2000a[ ). We investigated the magnitude of 
this effect in our data by altering the bandwidth and re- 
fitting the data. We adopted three different pass-bands 
in addition to the nominal 0.5-7.0 keV band, specifically 
0.7-7.0 keV, 0.5-2.0 keV and 0.4-7.0 keV. Although the 
response below ^0.5 keV is suspect, since most of the 
galaxies i n our sam ple were observed relatively early in 
the life of Chandra we expect that the degradation of the 



low-energy response, which is in part responsible for this 
uncertainty, should not be dramatic. Since the strongest 
O lines are in the range 0.6-0.8 keV and the S line is at 
~2.4 keV, we did not attempt to assess the impact of 
adopting the 0.7-7.0 keV or 0.5-2.0 keV bands, respec- 
tively, on the abundances of these species. 

First considering the 0.5-2.0 keV band, we found there 
are insufficient counts in the vicinity of the Fe L-shell 
clearly to disentangle the hot gas and point-source com- 
ponents in the faintest systems. In this case Zp c was typ- 



ically very poorly-constrained. In the brighter systems, 
the Fe L-shell is well-defined, which enables both tem- 
perature and metallicity to be determined reasonably, al- 
though degeneracy with the point-source component still 
persists, resulting in compromised error-bars and higher 
Zp c . Given the poor constraints wc did not include the 
magnitude of this effect in our systematic error tables. 

For the 0.4-7.0 keV band, we tended to find a slight 
increase in the measured Zp c , whereas in contrast, for the 
0.7-7.0 keV band, we found a significant reduction. Both 
of these results can be understood in terms of the impor- 
tance of constraining the continuum at energies below 
the Fe L-shell "hump" . If the continuum is not prop- 
erly estimated, the equivalent widths of the lines will be 
in error, resulting in an incorrect abundance determina- 
tion. Since the errors in the continuum measurement 
will impact all equivalent widths in the same sense, the 
abundance ratios are less affected by this uncertainty. 

5.5. Undetected sources 

We included a component to account for undetected 
point-sources in our sample. The adopted model has 
been shown to be a good approximation to the compos- 
ite spectrum of the detected point-so urces in a range of 
early- type galaxies ( Irwin et al. | |2003|) . It is by no means 
certain, however, that this model will provide a perfect 
description of the undetected X-ray point-sources in any 
galaxy. This especially may be the case since the spec- 
tra of fainter LMXB (which would be undetected) tend 
to be harder than those of brighter objects (Church 
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Balucinska-Church 2001 ). This difference may be exac- 
erbated if a substantial part of the bright X-ray binaries 
are high/ soft state black hole systems. 

To estimate the sensitivity of our results to this ef- 
fect, we adopted two tests. Firstly we allowed the tem- 
perature of the bremsstrahlung component to be fitted 
freely. In the brightest systems, this best-fitting temper- 
ature tended to be slightly lower than the canonical value 
(~2-3 keV). This may arise due to a slight inadequacy in 
our modelling of the hot gas; it tends to produce slightly 
higher best-fitting abundances, via a similar mechanism 
to the Fe bias. The abundance ratios were not signifi- 
cantly altered by this test. In the faintest systems, the 
temperature tended, if anything, to increase, which pro- 
duces a slight systematic reduction in the best- fitting Zp c . 
As a second test, we simply omitted the point-source 
component entirely. In the brightest systems, where the 
hot gas overwhelms the point-source contribution, the 
impact is negligible. In the fainter systems the omission 
of the point source component leads to a dramatically 
poorer fit and substantially lower abundances. This is 
easily understood as the continuum level is artificially 
raised to fit the high-energy residuals, systematically re- 
ducing the line equivalent widths. 

5.6. Hydrogen column- density 

Another source of systematic uncertainty in our anal- 
ysis is an error in our adopted hydrogen column to the 
galaxy. We experimented with allowing the Nn to fit 
freely for each galaxy. We found a strong anti-correlation 
between Zp c and iVu, since the abundance measurement 
is sensitive to the continuum at low energies. If Nh is 
under-estimated, the continuum level below the Fe L- 
shell will also be under-estimated, artificially raising Zp c . 
We found that the abundance ratios were relatively un- 
affected by this uncertainty. 

Typically the best-fitting Nb. was in good agreement 
with our adopted Galactic value. In a few cases the mea- 
sured value was slightly higher, although we do not be- 
lieve this indicates substantial intrinsic cold absorption 
in these galaxies, and is probably more representative of 
slight systematic uncertainties in the modelling and the 
responses at the lowest energies. 

5.7. Fe bias 

Although we have endeavoured to remove the Fe bias 
wherever possible, there exists the possibility that the 
quality of the data in some cases prevents us from de- 
tecting the presence of multi-temperature hot gas. We 
therefore include in our error-budget calculation a com- 
ponent to account for the Fe bias, if appropriate. In 
Table for those galaxies which require only one hot 
gas component, we show the impact of adopting a two- 
temperature model. For those galaxies requiring a two- 
temperature model, we show the impact of fitting a single 
hot gas component model. In order to make the sys- 
tematic error-bars as meaningful as possible in the total 
error-budget shown in Table [2J, however, we only include 
the Fe bias term for systems requiring just one hot gas 
model. 

6. DISCUSSION 
6.1. Near-solar Fe abundances 



In our sample of early- type galaxies, we did not find 
any convincing evidence for the very sub-solar Zp c histor- 
ically reported for such systems. For the brighter galax- 
ies, we found clear evidence that the ISM abundances are 
~solar. It is interesting to note that NGC 507 appears 
to be something of an "outlier" when compared to the 
other high-Lx systems, since its Zp e seems to be some- 
what lower than is typically seen in such objects. How- 
ever, at high Lx, the globally-averaged Zp c must start to 
fall in order ultimately to match Zp c ~0.4 typically seen 
in clusters, and since this is the highest-Lx object in our 
sample, its lower Zp e might be expected. 

Although the data for the fainter galaxies are typically 
of too poor quality to allow us to exclude very sub-solar 
abundances, they are generally consistent with Zp c ^l. 
In fact, if we consider all those systems for which only a 
single hot gas component was required in a single aper- 
ture (which are typically have poor er S/ N), excepting the 
possible outlier (NGC 1553; see § |6.2|) , which will bias 
the measurement due to its small error-bars, we obtain a 
mean Zp e =0.58tQ22- To investigate any possible correla- 
tion between Zp c and Lx or Lb we applied several statis- 
tical tests — Pearson's linear correlation test, Spearman's 
rank-order correlation test and Kendal's r test — to the 
data. In Fig | we show Zp c plotted against Lx and Lb- 
Since it is possible that the S/N of the faintest galaxies 
masks the presence of the Fe bias, in addition to test- 
ing the correlation with the whole dataset we separately 
considered a subset of the galaxies, in which spatially- 
resolved spectroscopy or two hot gas models were re- 
quired, thereby taking account of any possible Fe bias. 
Considering all the galaxies, Pearson's linear correlation 
test did not reveal a correlation between Zp e and either 
Lx or Lb (po, the probability of no correlation, being 32% 
and 63%, respectively). With the non-parametric tests 
we found evidence of a correlation with Lx, and marginal 
evidence of a correlation with Lb (e.g. for Spearman's 
test po=0.5% and 6%, respectively). Considering only 
the "reliable" subset of galaxies, however, all the tests 
failed to detect a correlation (p > 20%), suggesting that 
the correlations in the entire data-set may be artefacts 
related to the Fe bias, although this will need to be in- 
vestigated with higher-quality data. We conclude there 
is no convincing evidence of a correlation between Zp c 
and Lx or Lb- 

The principal reason for the under-estimate of Zp c re- 
ported in the literature appears to be overly simplistic 
spectr al modelling, as demonstrated by Buote fe Fabian 
( 1998 ). In many of the fainter systems X-ray point - 
sourccs c ontribute > 50% of the total X-ray flux. Using 



Chandrc we have been able to resolve a significant frac- 
tion of these point-sources, reducing the impact of this 
"contaminant" and allowing less biased abundance deter- 
mination. Simple fiducial models are frequently adopted 
to describe the spectral shape of the combined emission 
from undetected point-sources, but these have been by 
no means uniformly included in spectral-fitting in the 
literature. It is important to understand that these are 
determined empirically, and there is no a priori reason to 
believe them a perfect fit to the undetected source emis- 
sion in any given galaxy. This issue may have been, in 
part, responsible for the extraordinarily low abundances 
reported in three low-Ly gala xies observed with XMM by 
'Sullivan fe Ponman ( 2004 ) , since the flux from the hot 
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Fig. 5. — Best-fitting Zp c and 1-cr error-bars are shown as a function of Lx (left panel) and Lb (right panel). Data-points marked with 
a circles (black) represent those galaxies for which only a single aperture could be used, and only a single hot gas component was required 
by the data. The remaining galaxies are shown as triangles (blue). 



gas is most likely overwhelmed completely by the point- 
sources, making accurate abundance determinations very 
sensitiv e to this m odelling;. 

With Chandra\ we have also been able to resolve the 
spatially-varying temperature structure in a number of 
the brighter galaxies. The excellent agreement between 
our measure d abundances and prev ious single-aperture 
ASCA work (Buote & Fabian 1998), where the samples 
overlap, confirms that the temperature gradient is suffi- 
cient to explain the multiple hot gas components required 
in the (large) ASCA ex traction apert ures. A similar con- 



clusion was reached by Buote ( 1993) for a small number 
of bright galaxies, using Rosat, which we can now confirm 
with the su perior spe ctroscopic data, and finer spatial 
resolution of Chandra, in a wider range of galaxies. Re- 
solving the spatially- varying temperature structure mit- 
igates the "Fe bias", in which the addition of multiple 
components with different temperatures tends to sup- 
press the line equivalent widths, giving artificially lower 
abundances. A recent dramatic example o f this effect 



has a lso been seen in t he Antennae galaxies ( Baldi et al 
[2005 ), where Chandra\ has revealed complex ISM struc- 



ture in which Zp c is typically > 1, in contrast to previous 
(global) ASCA measurements of Zp c ~0.1. 

In general, we did not require multiple hot gas com- 
ponents in each annulus, indicating that the ISM is 
not multi- phase in these regions. Although, in most 
cases, the presence of multiple hot gas components in 
the ASCA apertures can be understood in terms of a 
temperat ure gradient which we are able to resolve with 



Chandra] , there still remains the intriguing presence of 
multiple hot gas components in single annuli for several 
of the galaxies. In most cases this is only required in 
the innermost bin, suggesting spatially-varying temper- 
ature structure at scales which were too small to ac- 
cumulate useful spectra. In NGC 1399 and NGC 5846, 
however, there was some evidence that multiple temper- 
atures may be needed over a range of annuli. In the 
case of NGC 5846, this simply seems t o reflec t the r ather 
disturbed morphology of the galaxy. Buote ( 2002| ) first 
commented on the "limited multiphase gas" in NGC 1399 
based on XMM data, and a similar feature is also seen 
in several groups (NGC 5044: |Buote et al.| |2003a| , M 87: 



pastaldello fc Molcndi| |2002; , RGH80: |Xuc et al.||2004| ). 
This two-temperature structure may be related to the 
heating; of the gas an d the quenching of the cooling flow 
flMathews et al.||2004j ). 

It is worth noting that the shapes of the abundance 
profiles we measured are consistent with ob servations of 
the ce ntre of the bright group NGC 5044 (|Buote et al 



2003b), although we were only usually able to obtain 
constraints for the most massive systems, e.g. NGC 1399. 
Nonetheless, these abundance gradients are similar to 
what are predi c ted b y the "circ ulation flow" model of 
Mathews et aL] fl2004j) (sec § O). 



6.2. Stellar versus ISM abundances 



Arimoto et al. ( 1997 ) pointed out that the historically 
low ISM abundances for early-type galaxies were in stark 
disagreement with the abundances determined for the 
stellar content. This is difficult to interpret in terms of 
simple chemical enrichment models in which the gas is 
enriched by stellar mass-loss and supernova ejecta. In 
such a scenario, we might expe ct t he ISM abundances to 
exceed those in the stars (see § 6.4 for more discussion of 
these models). In several galaxies we found best-fitting 
abundances which were moderately sub-solar (~0.5), and 
in one case (NGC 1553) we did, in fact, obtain a best- 
fitting Zp which was ~0.1 (although adopting a two- 
temperature hot gas model, which may be more realis- 
tic in this disturbed system, the data were not inconsis- 
tent with approximately the same ISM and stellar abun- 
dances). To assess quantitatively whether there is still 
any evidence of the discrepancy between stellar and ISM 
abundances we searched the literature for optical line in- 
dex measurements, from which we were able to estimate 
the global Fe abundances in the stars. We discuss this 
analysis in detail in Appendix |A|. Table A5 



summarizes 



the properties of the stellar population of each galaxy. 

We show a comparison of the stellar and ISM Zp e in 
Fig [| The majority of the data-points cluster close to 
the ZFc(gas)=ZFo(stars) line. Clearly a number of points 
lie significantly above this line, suggesting that these 
galaxies have experienced significant ISM enrichment. A 
few points lie slightly below the line, but are approxi- 
mately consistent with it within their error-bars. Only 
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Fig. 6. — A comparison of Zp c in the hot gas versus estimated mean stellar Fc abundance. We show separately the data-points for which 
spatially-resolved spectroscopy or multiple hot gas components were required (triangles) and those for which only a single temperature, in 
a single aperture was needed (circles) . When no upper limit could be obtained for any galaxy, we show only the lower limit (arrows) . We 
show the line indicating equal metallicity in both components. For clarity, any data-points lying below Zp e = 0.1 are artificially raised to 
this level and all error-bars shown are 1-cr. 



NGC 1553 lies significantly below the line and this galaxy 
has a disturbed morphology, so that there is considerable 
uncertainty in using only one hot gas component to fit 
the spectrum. Alternatively, the ISM abundance for this 
object may in fact be low, in which case the galaxy must 
have acquired a substantial amount of uncnrichcd, pri- 



With our Chandra data we have been able to obtain, 



more ial gas. Why such a process should seem to operate 



so strongly only in this galaxy, and whether it is related 
to the disturbed morphology are not entirely clear. 

If we consider those galaxies for which an upper-limit 
on the abundance could not be obtained, we find that 
the best-fitting values (not shown in Fig [|) of 4 of these 
systems lie above and 8 below the line. This is what we 
would expect if the distribution of Zp e in these systems 
is comparable to the "detected" galaxies. We therefore 
conclude that these galaxies probably have a similar Fe 
enrichment history, although this remains to be investi- 
gated with hig her-qualit y data. 

Based on the Chandra data, therefore, we found no ev- 



idence that the ISM abundances of the early-type galax- 
ies in our sample (with one exception) are significantly 
lower than the mean metallicity of the stars. Further- 
more, many of these galaxies may have ISM abundances 
slightly higher than the stellar component. 



6.3. Supernova enrichment 



in many cases for the first time, interesting constraints 
on a-to-Fe abundance ratios. A number of competing ef- 
fects can enrich or dilute the ISM metal content in early- 
type galaxies, including supernovae ejecta, st ellar winds 
and large-scale g alactic gas inflow/ outflow (Matteucci 
fe Grcggkj 1986f ), The a-elements and Fe are predomi- 



nantly processed through supernovae, so stellar mass-loss 
acts only as a "slow-release reservoir" for the elements; 
if the stars were all formed very quickly, we would ex- 
pect the metal yields of the mass-losing stars to resemble 
SNII. Because these element are primarily produced in 
supernovae it is common-practice to define, as a fiducial 
diagnostic, the fraction of the Fe content of a galaxy pro- 
duced in SNIa, f. Since Fe is more profusely synthesized 
in SNIa than SNII, typically this quantity is measured 
by matching linear combinations the SNe yields for a 
complete stellar population (integrated over an assumed 
IMF; typically that of S alpeter) to the abundance ra tios 
with respect to Fe (e.g. |Gastaldello fc Molcndi| |2002| ) . 



It is important to stress that, on account of the finite 
timescale over which the gas is processed through stars 
the SNIa fraction measured in the hot gas, fg| s and that 

determined from the stars, f* a , are not identical, and nor 
are they entirely commensurate with the total fracti on of 



Fe synthesized in SNIa during the galaxy history (Mat 
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TABLE 4 

Supernova enrichment fractions 



Galaxy 


Elements 


All elements 








Omitting 


o, s 










yield 


X 2 /dof 


pla 

gas 


r la 


;e 


yield 


X 2 /dof 


rla 
gas 


rla 

i rangc 


High-Lx galaxies 


IC4296 
NGC 507 


0,Mg,Si 

0,Ne,Mg,Si,S,Ni 


WDD2 
W7 


0.035/2 
16/5 


0.60 ±0.16 
0.67 ±0.05 


0.28- 
0.46- 


-1.0 
-1.0 


WDD2 
W7 


0.006/1 
0.95/3 


0.60 ±0.16 
0.70 ± 0.05 


0.28-0.63 
0.46-0.98 


NGC 741 
NGC 1132 


0,Mg,Si,Ni 
0,Mg,Si,Ni 


W7 
WDD1 


8.5/3.0 
11/3 


0.88 ± 0.05 
0.87 ±0.07 


0.75- 
0.65- 


-0.91 
-0.87 


W7 
W7 


1.6/2 
3.9/2 


0.70 ±0.13 
0.66 ±0.12 


0.45-0.94 
0.43-0.77 


NGC 1399 


0,Ne,Mg,Si,S,Ni 


WDD2 


180/5 


0.87 ±0.01 


0.70- 


-0.91 


W7 


3.4/3 


0.74 ± 0.02 


0.49-1.0 


NGC 1600 
NGC 4472 


0,Mg,Si 

0,Ne,Mg,Si,S,Ni 


W7 
WDD1 


0.18/2 
140/5 


0.68t°;« 
0.81 ± 0.02 


0.37- 
0.53- 


-0.96 
-0.81 


W7 
W7 


0.16/1 
11/3 


fi7+0- 16 
U - D '-0.21 

0.60 ± 0.03 


0.37-0.75 
0.31-0.66 


NGC 5846 
NGC 7619 


0,Ne,Mg,Si,S,Ni 
0,Mg,Si,Ni 


WDD2 
WDD2 


58/5 
24/3 


0.88 ±0.02 
0.78 ±0.06 


0.71- 
0.49- 


-0.92 
-0.88 


W7 
W7 


2.0/3 
4.4/2 


0.73 ± 0.04 
0.63 ±0.07 


0.49-0.87 
0.27-0.73 


NGC 7626 


Mg,Si 










W7 


0.004/1 


0.53 ±0.29 


0.17-0.59 


JVloderate-Lx 


galaxies 




















NGC 720 1 


0,Ne,Mg 


WDD1 


21/2 


0.85 ±0.06 


0.63- 


-0.85 


WDD1 


2.4/1 


0.58 ± 0.12 


0.02-0.58 


NGC 1332 1 


0,Ne,Mg,Si 


WDD2 


37/3 


0.72 ±0.05 


0.39- 


-0.75 


W7 


2.7/2 


0.55 ± 0.07 


0.16-0.56 


NGC 1387 
NGC 1407 


0,Mg,Si 
0,Mg,Si,S,Ni 


W7 
WDD1 


0.003/2 
39/4 


04+0.O6 
0.84 ± 0.04 


0.82- 
0.61- 


-0.94 
-0.84 


W7 
W7 


0.002/1 
1.8/2 


n 04+0.O6 
0.58 ±0.08 


0.82-0.94 
0.24-0.62 


NGC 1553 


O.Ne 


WDD2 


13/1 


0.71 ±0.09 


0.33- 


-0.71 


WDD1 


0/0 


0.42 ± 0.16 


0.0-0.43 


NGC 1700 


0,Ne,Mg 


WDD2 


3.2/2 


0.84 ±0.08 


0.65- 


-0.84 


WDD1 


2.3 x 10- 7 /l 


0.76 ±0.11 


0.49-0.76 


NGC 3607 
NGC 3923 


Mg 

0,Ne,Mg,Si 


WDD1 


39.1/3 


0.89 ± 0.05 


0.71- 


-0.97 


W7 
WDD1 


0/0 
5.4/2 


0.76 ± 0.24 
0.56 ±0.11 


0.33-0.77 
0-0.56 


NGC 4365 


0,Mg 


WDD1 


0.48/1 


n 77 +o.i9 
u.( <_ 27 


0.50- 


-0.77 


WDD2 


0/0 


0.67 ±0.34 


0.076-0.67 


NGC 4552 


0,Ne,Mg,Si,Ni 


WDD1 


76/4 


0.89 + 0.02 


0.73- 


-0.92 


W7 


3.1/3 


0.73 ± 0.04 


0.42-0.74 



Note. — The inferred SNIa enrichment fraction, fg% s , for the ISM of the sample galaxies. The preferred SNIa metal yield model is 
listed for each galaxy ("yield"), whereas we adopt the N97 metal yields for SNII. The range of best-fitting f ^ values, f T ^ n „ e , f° r different 
combinations of SNII and SNIa yields (sec text) are also shown. We show the results derived from all the measured abundance ratios ("All 
elements") and those which omit, if possible, O and S, both of which seemed aberrant (sec text). 



teucci & Chiappini 2005) 



However fg| s and fl a 



can be 



calibrated ag ainst realistic gasdynamical e nrichment sim- 
ulations (e.g. Brighcnti & Mathews 1999| ), making them 
powerful measures of the enrichment history of a galaxy; 
in any case, they provide a useful yard-stick with which 
to compare different galaxies. Furthermore the relative 
abundances of elements which are predominantly prod- 
ucts of one type of supernova are fixed by the SNe yields, 
since there is unlikely to be significant segregation of met- 
als produced together. This enables measurement of the 
abundance ratios to provide direct constraints on the SNe 
yields and, in principle, the IMF. Given the remarkable 
degree of scatter in the pr edicted metal yield s between 



different explosion models ( |Gibson et al 



1997, hereafter 



G97)i , such constraints are of particular importance 



In I order to determine values of ill 



. -gg fo r each galax y- 

we p erfo r med fi tting analogous to that of Gastaldello & 



Molcndi ( 2002 ). We experimented with different SNe 
metal yields, which are highly sensitive to the supernova 
modelling. We adopted th e SNIa W7, WDD1 or WDD2 
yields from Nomoto ct al. ( |l997b| ) and, for ease of com- 
pariso n with other authors, the SNII yields from Nomoto 
et al. (1997a, hereafter N97). We also experimented with 
a range of SNII yields taken from G97, in order to exam- 
ine the sensitivity of the results to the adopted model. In 
practice the G97 yields are less useful than N97 since they 
do not all include the same range of species. It is worth 
noting that all of these models assume a Salpeter IMF, 
although the IMF is not expected to affect the abundance 
ratios as significantly as the overall abundances (G97). 
We note that fitting this model to the Zo/Zp e , Zno/Zfc, 
ZMg/ZFc, Zsi/Zpc, Zs/Zp and ZNi/Zp c ratios observed 



in the Sun (Asplund ct al. 2004) yields the well-known 
fj: a ~0.75 appropriate for the Solar neighbourhood whilst 
strongly favouring the WDD2 SNIa yields. 

Initially we fitted our model to all of the available a /Fe 
abundance ratios, adopting the SNIa model which best 
fitted the data. There is remarkable agreement between 
the supernova fractions measured in each galaxy, with 
fga S ~ ranging from 0.7-0.9, although it fell as low as ^0.3 
if the G97 yields are used. The results are summarized 
in Table ^. Such values are in excellent agreement with 
recent observations of the centres of groups and clusters 



For example f g | s -0.8-0.9 in M 87 ( |Gastaldello & Molendil 



2002), fgL-0.7-0.8 in NGC 5044 (Buote ct al. 2003b) and 



rla r 

gas 



0.8-0.9 in the centre of the cluster A 1795 (Ettori 



et al. 2002 ). We found a general trend that the delayed 



detonation SNIa yields were preferred over the deflagra- 
tion (W7) model, although this seems to be driven by 
the low Zo/ZMg ratios in these models, which could be 
a source of uncertainty (see below). 

In contrast to the enrichment of the gas, super-solar 
a/ Fe ratios are typically reported in the stars when sim- 
ple stellar population models are fitted (see Appendix^), 
which implies substantially more Type II enrichment m 
the stars. Assuming the a/Fe abundances obtained in 



Table A5 are driven primarily by the Mg/Fe ratio (and 
assuming that there is no stellar Mg/Fe gradient), we es- 
timate an average fjf~ 0.35, with individual values rang- 
ing from - 0.60 for NGC 5018 to - 0.1 for NGC 1399, 
adopting the W7 SNIa and N97 SNII yields. This is 
a well-known result (e.g. Worthcy 1998| ) which is typi- 
cally explained in terms of the short timescale of star- 
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formation in early-type galaxies, or variations in the 
IMF. The stark discrepancy between fg| s and il a 



can 



readily be explained only if substantial amounts of Fe 
are added to the ISM through SNIa occurring after the 
bulk of the stars have formed. Such a large input of Fe 
into the ISM would tend to raise Zp e to ~2-3 times the 
mean stellar value. The error-bars on our measured Zp c 
in the gas are sufficiently large that many of the galax- 
ies in the sample could have ZpoConsistent with such an 
elevated value; in particular we note that NGC 1332 and 
NGC 4472 have best-fit abundances close to that required 
by this argument. However, several galaxies (most no- 
tably NGC 1399 and NGC 720) seem to have Z Fe about 
a factor ~2 too low, implying dilution of the ISM by a 
source of primordial (i.e. un-enriched) gas in these sys- 
tems. 

Despite consistency with f g | s determined in other sys- 
tems, it is clear that the fits are typically rather poor, 
suggesting a disagreement between the standard SNe 
yields and the ISM abundance pattern. This primarily 
arises because of discrepant a-element abundance ratios, 
even taking into account the large scatter in the SNII 
yields. The chiefly discrepant species is O, which ap- 
pears to be significantly under-abundant in comparison 
to other a-elements, such as Mg (for which we have the 
most measurements) . The measurement of reliable Zo in 
the X-ray band is challenging, especially when attempt- 
ing imaging spectroscopy at energies where the CCD re- 
sponse may not be perfectly well understood (§ o). At 
plasma temperatures ~0. 5-1.0 keV, appropriate tor our 
sample, the principal O features tend to be blended some- 
what with the Fe L-shell "hump" which exacerbates the 
difficulty in reliably determining Zo- However, unexpect- 
edly low values of Zo/ZM g (~0.8-1.3, given our default 
solar abundances) have been widely reported in the liter- 
ature, in contrast to Zo/ZM g ~l-4-2.8 predicted by SNII 
models (G97). This discrepancy has been found in both 
imaging and grating spectroscopy data from a variety of 
different instruments observing a wi de range of systems 



such as the bright g roup N GC 5044 (Buote et al. 2003b 
Tamura et af1|2003), M87 flGastaldeflo k Molend]| 



gakelliou et alj p002|) , the massive elliptical NG C4636 
( Ku et al. 2002| ), the starbursting galaxy M82 flTsuru 



et al _ 1997 ) and th e centres of some galaxy clusters ( [Pe- 
terson et al. |2005). There is even some evidence of thTs 



effect in the Milky Way ISM flUeda et alj p005| ). 

It is not certain whether the observed O deficiency rep- 
resents a problem with the SN yields, which are highly 
sensitive to the (rather uncertain) SNe modelling (G97), 
or if it is indicative of an additional source of metal en- 
richment within the galaxy. It is interesting to note that 
some SNR have been observed with apparent abundance 
anomalies (for example N49B in which the ejecta appear 
signifi cantly enriched in Mg without accompanying Ne 



or O: Park et al 



20031), which might hin t at a prob- 
l em w ith the nuclear physics. In contrast, Loewcnstcin| 
( pOOl ) suggested there may be an alternative source of 



metal enrichment to explain the low Zo/Zpo in clusters. 
He proposed that significant enrichment in the early uni- 
verse may arise from population III hypernovae, since the 
O-burning region is significantly expanded in an hyper- 
nova. Hypernova m etal yields have been c o nrput ed by 



unfortunately they did not exhibit the very low Zo/ZMg 
(nor Zo/Znc) which are observed in our data, although 
the calculations for Mg and Ne production are sensitive 
to the model ass umptions. Given the publi shed Pop III 
hypernova yields Matteucci & Pipinc (2005) pointed out 



they would not have a significant impact upon the abun- 
dance pattern of an elliptical galaxy. 

Given these reservations, we did not attempt to in- 
clude a putative hypernova component in our modelling. 
If, instead, the discrepancy between the observed and 
predicted abundance of O reflects a problem in the com- 
puted yields, notwithstanding their success in reproduc- 
ing the Solar abundance pattern, an alternative approach 
is simply to omit this element from the computation of 
fg| s . This tended to result in slightly lower estimated val- 
ues of fg|g, but substantially improved \ 2 - Interestingly 
we also found that Zs/Zp e , where it could be measured, 
was slightly under-pred icte d by the model (which may 
be a systematic error; § 5T) and so we also omitted this 
element, although it had little impact on the best-fitting 
f g | s . The fit results are shown in Table |]. We found 
a mean fi| s =0.66 ± 0.11, where the quoted error is the 
standard deviation in the mea sure ments . Using the same 
statistical tests described in 



lations between fi* and Lx 

gas ^ 



Umeda ct alj ( |2002j ) and |Heger fe Woosley| (|2002J) , but 



6.1 we searched for corre- 
and Zp e , but found no 
evidence of any such correlations (po > 40% in all cases) . 
In this case, we found that the W7 (deflagration model) 
abundance yields are preferred over the delayed detona- 
tion models, although in most cases (with one or two 
notable exceptions, such as NGC 1399) the data do not 
strongly discriminate between them. It is striking that 
f g | s in these galaxies is remarkably similar to that in the 
Solar neighbourhood, implying a similar enrichment his- 
tory for the hot ISM of early-type galaxies and the cold 
ISM of spirals. 

An intriguing piece of evidence which further compli- 
cates this picture of metal enrichment is the observation 
that [O/Mg] in Galactic bulge st ars appears to be a de - 
creasing function of metallicity (Fulbright et al. 2004| ) . 
There is a reduction in [O/Mg] by as much as ~0.5 dex 
with respect to the SN II yield level by Zp c ~l, providing 
evidence that O may be under-abundant in old stellar 
populations, not just the ISM. Moreover such a trend is 
very difficult to reconcile with a picture where only SNII 
and SNIa enrich the gas. Nonetheless whether a putative 
Pop III hypernova population could resolve this issue is 
unclear. 



6.4. Enrichment models 

We now discuss our results in the light of various mod- 
els for the chemical enrichment of early- type galaxies. 
In traditional "galactic wind models" the galaxies form 
monolithically and passively evolve, during which time 
the dynamics and luminosity of the gas is strongly af- 
fected by supernova feedback, which can drive pe riods of 
large scale galactic outflow (e.g. Giottietal. 1991). Given 
the large amount of metals injected into the ISM by 
the supernovae, these models tend to predict Zp e ~a few 
times solar, (and, certainly, dramatically in excess of the 
abundance of the stars), which is clearly in disagreement 
with the majority of galaxies in the sa mple, in many case s 
by as much as an order of magnitude. Pipino ct al. ( 2005 ) 
allowed the accretion of additional, unenriched gas, but 
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were unable to mitigate this discrepancy while reproduc- 
ing the oth er optical and X-ray propert ies of the gas. 
In contrast, Brighcnti fe Mathcwsj (119991 ) constructed a 
cooling- flow gasdynamical model for NGC4472, incor- 
porating the prolonged infall of gas. These authors were 
able to produce ~0.5 solar abundances for the galaxy (the 
precise abundances being rather sensitive to the assumed 
supernovae rates), although the density profile was too 
peaked on account of excessive cooling in the core. Fur- 
thermore, the enrichment history assumed in this model 
makes it difficult to understand the large amount of met- 
als in clusters in terms of the contribution from elliptical 
galaxies 



An alternative approach was adopted by Kawata & 
Gibs^ d2003|) , who followed the evolution of a ~ 1O 13 M 



elliptical galaxy in an hierarchical structure formation 
simulation. In their model, however, most of the Fe pro- 
duced by the stars is locked into large amounts of cool 
(~ 10 4 K) gas and excessive on-going star formation, 
at variance with observations. In this case the hot phase 
ISM was found to have Zp e < 0.3 and, critically Zp e (stars) 
was ~an order of magnitude greater than that of the hot 
gas. The results are clearly at variance with our Chandrc\ 
measurements, in the opposite sense to the galactic wind 
models. 

One key ingredient not typically incorporated into en- 
richment models is the role of AGN feedback. There is 
some evidence of disturbances in some "normal" early- 
type galaxies which may point to interaction between the 
ubiq uitous supermassive central black-hole and the ISM 



(e.g. Blanton ct al. 2001; Jones ct al. 2002| ). Numerical 



simulations have suggested that ISM- AGN interactions 
can have a dra matic effect on the dynamics and energet- 
ics of the ISM (|Di Matteo et al.||2005|) . |Kawata fe Gib" 
son[~( 2005| ) incorporated AGN feedback into their model, 



finding that this effectively evacuated the (cool) enriched 
gas from the galaxy, so that the discrepancy between 
stellar and ISM abundances of their model persists. We 
note that these authors comment on the good agreement 
between their model an d Zp P for the galaxy NG C 3923, 
based on the results of 



Matsushita ct al.| (|2000| ). O 



analysis of the Chandrc data indicates a significantly 
higher abundance for this galaxy, underlining the dis- 
crepanc y between our r e sults and the predictions of this 
model. Mathews et al. ( 2004 ), in contrast, showed that 
allowing a small amount of continuous heating from a 
central AGN may drive large-scale "circulation flows", 
which efficiently mixes gas through the galaxy. This 
model could repr oduce th e abundance gradients found 
in the XMM and \Chandrc\ data of the group NGC 5044. 
Crucially it predicted the correct Zp c , which is coinci- 
dentally similar to that observed in most galaxies in our 
sample. Whether a variant of this model can successfully 
describe the lower-mass systems studied in this present 
work remains to be seen, however. 

7. SUMMARY 

We have surveyed the metal abundances of the hot 
ph ase ISM w ithin a sample of 28 early-type galaxies, us- 
ing 



Chandra. In summary: 



there convincing evidence that Zp e correlates with 
the galaxy luminosity. Considering just our lowest- 
S/N galaxies, we found a mean Zf o =0.58+q;22. 
The discrepancy with past results appears to be 
due to imperfect modelling of t he (unde tected) 
point-source populations (of which Chandrc has en- 
abled a substantial fraction to be detected) and 
the fitting of single-temperature models to multi- 
temperature gas (the Fe bias). 



2. We confirmed the conclusions of Buote & Fabian 



(1998) that multiple components are often required 
to fit the spectra of early-type galaxies. In anum 



ber of cases the excellent spatial resolution of Chan 



1. We did not find the historically very sub-solar 
abundances reported for early-type galaxies, nor is 



drc enabled us to resolve temperature gradients ly- 
ing within ASCA or Rosat spectral extraction aper- 
tures. In a few cases the fit required multiple hot 
gas components in a given annulus. In all cases a 
component was required to account for undetected 
point-sources. 

3. Comparing the emission- weighted ISM and stellar 
Fe abundances we found excellent agreement. Ex- 
cepting one possible outlier (NGC 1553) there is 
no evidence that the stellar abundances are signifi- 
cantly higher than those of the gas. In fact the data 
allow the ISM Zp e to be slightly higher on average. 

4. The data are inconsistent with both monolithic 
"galactic wind models" , which ov er-predict Zp and 
the re cent hierarchical models of Kawata fe Gibson 
( |2005 ), which under- predict Z p e. Th e "circulation 
flow" model of Mathews et al. ( [2004 ) is able to ex- 
plain Zp c of about the correct magnitude, albeit for 
a system more massive than most considered here. 

5. We found that the SNIa Fe enrichment fraction in 
the hot gas, fg| s , is 66 ± 11% and shows no obvious 
trend with luminosity. This value is remarkably 
close to the 75% enrichment fraction for the solar 
neighbourhood, indicating similarity in the enrich- 
ment histories for the early-type and spiral galax- 
ies. It is also similar to that observed in the cen- 
tres of clusters and groups, indicating homology of 
enrichment from cluster scales to intermediate-Lx 
galaxies. 

6. The ZMg/ZFe ratio in the gas is significantly lower 
than that in the stars, indicating that significant 
SNIa must enrich the gas, in addition to stellar 
mass-loss. In at least two cases (NGC 1399 and 
NGC 720) this would seem to require the further 
accretion of unenriched, "primordial" material to 
prevent Zp e being significantly higher than ob- 
served. 

7. There is increasing evidence that O is over- 
predicted by the SNII yields, which may imply 
problems in the supernova calculations. Alter- 
natively, the measurement of varying [O/Mg] in 
Galactic bulge stars may point to a source of a- 
element enrichment in addition to SNIa and SNII. 
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APPENDIX 



STELLAR ABUNDANCE ESTIMATES 

In o rder to estimate the mean Fe abundance in the stars, we searched the literature to obtain reliable Lick/IDS 
index (Faber et al. 1989) measurements for the centre of each galaxy. In order to break the age-metallicity degeneracy, 
we converted the line indices into [Fe/H] by linearly interpolating the simple stellar population (SSP) model results of 
Thomas et al.|( ^003[ ). T hese models take in to account non-solar a-element to Fe ratios, which are typically reported in 



early-type galaxies (e.g. [Trager et al.||2000| ). Where possible, we fitted the interpolated model with a x minimization 
procedure to the central H/^, Mgb, Fe5270 and Fe5335 indices. In several cases some or all of these indices were not 
available in the literature, in which case we either adopted appropriate alternative indices or omitted them. In a 
number of cases one or more of the indices appeared to be discrepant with the others, which may be a consequence 
of systematic measurement errors (or, possibl y, t he breakdown of the SSP model assumptions). In these cases, we 
omitted the most discrepant index. In Table A5 we list the indices used for each galaxy. For NGC 1387, we could 
not find usefu l line index measurements in the literature, and so we omitted it from the comparison with the X-ray 
results (§ |6.2|) . For several galaxies, it was necessary to fix the age or [a/Fe] in order to obtain interesting abundance 
constraints. In these cases, we adopted an age of 12 Gyr and [a/Fe]=0.23, which are roughly consistent with the 
average valu es found in the r emain ing galaxies. The best-fitting parameter values, converted to the solar abundances 



standard of Asplund et al 



(2004), and the reference from which the line strength measurements were derived are 
shown in Table A5. 

It is interesting to compare our results with tho se of other recent ca talogues of ages and metal abundances in 



early-type galaxies which overlap with our sample. Thomas et al. ( 2005 ) measured ages and central abundances for 
a sample of galaxies, of which 14 overlapped our sample, using the SSP models adopted in this present work. In 
general we found good agreement with our results; the only two galaxies for w hich Zp r was substantially dis crepant 
were NGC 1 600 and NGC 7619, for which our abundances were ~0.1 dex higher. Terlevich & Forbes (2002) and Trager 



Est al.| ( 2000) both est i mate d (central) metal abundances, ages and a/Fe ratios using the older SSP models of Worthey 
( |1994D !~ [rrager et alj (20001) made c orrections for non-solar a/Fe ratios so that [Fe/H] could explicitly be obtained. In 
contrast 
([Z/H] 



Terlevich & Forbes ( 2002| ) adopted the [MgFe] lick index, which is more sensitive to the overall abundance 
[Fe/HJ + 0.94[a/Fc[ for the models we adopt). In general our results agree reasonably well (±2cr) with these 
previous measurements where the data overlap. There was som e scatter, however; for example we found that our 
central [Fe/H] were typically ~0.1 dex higher than th ose found by | Trager et all ( 200C ). We found the most significant 
discrepancies (> 3<r) between our results and those of Terlevich & Forbes "( 2002[ ) for the galaxies NGC 1553, N GC 3115 
and NGC 36 07, for all of which those authors reported poor quality data. Given the good agreement with Thomas 



|at al. ( 2005 ), the differences are likely a consequence of the different SSP models used. 



To estimate the total Fe content of the stars it is necessary to take into account any abundance gradients. T here 



is a substantial body of literature on line index gradients in early-type galaxies (see Kobayashi & Arimotc 1999| , and 
references therein) . Probably the most extensive collection of a bundance gradient measurements in early-type galaxies 
overlapping our sample is that of Kobayashi & Arimoto ( 1999|) who used linear approximations for the line-strength- 
metallicity relations at fixed galaxy ages and treated each line index separately. Unfortunately, their results were very 
sensitive to the line index adopted, and their method for relaxing the age constraint was not always successful for each 
galaxy (see their results for NGC 720) . A fully self-consistent assessment of stellar abundance g radients in al l of ou r 



galaxies is beyond the scope of this present work. However, for those galaxies in the sample of Trager et al 
we fitted simultaneously the data in both the r e /8 and r c /2 apertures (where r is the effective radius of the galaxy), 
assuming that the age and a/Fe ratios are th e same in each aper ture. This allowed us to estimate the magnitude of 
the stellar abundance gradient. We note that [Trager et al. ( 2000 ) allowed for both age and a/Fe ratio gradients, but 
we typically did not find this necessary to obtain good fits. In any case, if the overall stellar population comprises 
a mixture of different-aged components, the application of single burst SSP models is probably overl y simp listic. 



To estimate the global emission- weighted Fe abundance, we adopted the prescription of [Arimoto et al.| ( 1997 ), who 
assumed Zp c oc r _c for the stars, where r is the distan ce from the gala ctic centre. Assuming a circular aperture, which 
is appropriate for these line index measurements (see Gonzalez 1993), and a de Vaucouleurs optical brightness profile, 
we can express the emission weighted abundance within aperture r as: 

7 - M 7 - , N P(8-4c,7.67(r/r e ) 1 / 4 ) 
Z Fe (r) = Z Fe {oo) p(8;767(r/re)1/4) 

where P(a,x) is the incomplete gamma function J* dt cxp(— t)t a_1 /T(a). Two abundance measurements in different 
apertures are therefore sufficient to constrain both unknowns, c and ZFe(oo), of which the latter is the quantity of 
interest for comparison with the X-ray results. Following convention, we henceforth denot e log 10 (Zp e (r e /8)) = [Fe/H]o 



al 



and log 10 (ZF O (oo)) =< [Fe/H] >. For the 10 of our galaxies in the Trager et 
weighted mean of the difference between < [Fe/H] > and [Fe/H]n was 
in good agreement with the average ~0.3 found by Kobayashi fe Arimotc ( 1999 ) 



(200C) sample, we found that the 
0.27 -jO.06. This corresponds to c=0.24, 
Therefore, we used this factor 



to estimate the global abundances where we only had central line indices. The assumed power law relation for Zp c 
must break down at small radii, since it inflates unphysically, and so, for convenience, we assume that any "central" 
abundances inferred from the literature are the same as would be measured in an r e /8 aperture. In practice correcting 
from an aperture of size r e /10 to one of size r c /8 involves shifting the abundances only by 0.02 dex, which is considerably 
smaller than the statistical errors. 
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TABLE A5 
Stellar population parameters 



Galaxy 




indices 


rcf. 


[a/Fe] 


age 


[Fe/H] 


< [Fe/H] > 
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Mgb, 
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5 


0.20±0.1o 


12 


0.09±0. 22 


-0.13±0.23; 


INGC 1132 
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! 




3 


0.23 


12 


-0.25±0.11 


-0.47±0.13| 
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H/3, 


Mgb, 
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i 
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19+4 
12_ 3 
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-U.ll±U.lo!£ 


NGC 1407 


H/3, 


Mgb, 


Fe5270, Fe5335 


1 


0.33±0.02 


12±2 


0.04±0.06 


-0.18±0.08; 


NGC 1549 


H/3, 


Mgb, 


Fe5270, Fe5335 


1 


0.24±0.01 


5.1±0.7 


().12±0.03 


-0.10±0.07| 


NGC 1553 


H/3, 


Mgb, 


Fe5270, Fe5335 


1 


0.17±0.01 




0.23±0.04 


0.01±0.07t 


NGC 1600 


H/3, 


Mgb, 


Fe5335 


4 


0.24±0.02 


6.6±2.3 


0.30±0.07 


0.37±0.21 


NGC 1700 


H/3, 


Mgb, 


Fe5270, Fe5335 


4 


0.16±0.02 


n g + 0.5 

^•°-0.4 


0.31±0.08 


0.16±0.12 


NGC 3115 


h/3, 


Mgb, 


Fe5270, Fe5335 


5 


O.lliO.ll 


12 


0.25±0.16 


0.03±0.17t 


NGC 3585 


H/3, 


Mgb, 


Fe5270, Fe5335 


5 


0.20±0.12 


12 


0.09±0.17 


-0.13±0.18t 


NGC 3607 


H/3, 


Mgb, 


Fe5270, Fe5335 


5 


0.19±0.12 


15±4 


0.09±0.17 


-0.13±0.18j 


NGC 3608 


H/3, 


Mgb, 


Fe5270, Fe5335 


1 


0.19±0.03 


8.2±2.0 


0.16±0.07 


-0.17±0.16 


NGC 3923 


H/3, 


Mgb, 


Fe5270, Fe5335 


1 


0.34±0.03 


n q + 0.7 

°-°-0.3 


0.35±0.08 


0.13±0.10t 


NGC 4365 


H/3, 


Mgb, 


Fe5270, Fe5335 


5 


0.19±0.14 


12 


0.20±0.18 


-0.02i0.19j 


NGC 4472 


H/3, 


Mgb, 


Fe5270, Fe5335 


4 


0.25±0.03 


9.0±2.0 


0.13±0.06 


-0.01±0.20 


NGC 4494 


H/3, 


Mgb, 


Fe5270, Fe5335 


5 


0.16±0.13 


12 


-0.02±0.19 


-0.24±0.19t 


NGC 4552 


H/3, 


Mgb, 


Fe5335 


4 


0.24±0.02 


12±1.7 


0.17±0.04 


-0.05±0.10 


NGC 4621 


H/3, 


Mgb, 


Fe5335 


5 


0.28±0.14 


12 


0.19±0.17 


-0.03±0.18t 


NGC 5018 


H/3, 


Mgb, 


Fe5335 


5 


0.01±0.14 


z - u -0.4 


0.37±0.27 


0.15±0.28t 


NGC 5845 


H/3, 


Mgb, 


Fe5335 


5 


0.26±0.17 


12 


0.10±0.25 


-0.12±0.26t 


NGC 5846 


H/3, 


Mgb, 


Fe5335 


4 


0.22±0.02 


15±2.0 


0.05±0.04 


-0.18±0.16 


NGC 7619 


H/3, 


Mgb, 


Fe5335 


4 


0.18±0.02 


13.5±1.0 


0.21±0.04 


-0.12±0.15 


NGC 7626 


H/3, 


Mgb, 


Fe5335 


4 


0.24±0.02 


14±2 


0.08±0.05 


-0.25±0.18 



Note. — The mean abundances determined from fitting SSP models to the Lick indices reported 
in the literature. Those mean stellar abundances (< [Fe/H] >) marked | were estimated from 

otherwise the mean 
the parameter wat 
LEDA A— Tragei 



the central abundance ([Fe/H]o) assumin, 
abundance is inferred by extr. 



et al 



Table reference 
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